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The  techniques  of  pattern  recognition  are  proposed  as 
being  able  to  bring  together  any  analytical  problem  with 
an  existing  suitable  analytical  technique.   An  example  of 
the  analytical  capabilities  of  pattern  recognition  is  given 
via  the  descrimination  of  the  geographical  origins  of 
frozen  concentrated  orange  juice  by  looking  at  the  trace 
metal  content  of  five  selected  fingerprinting  elements 
(barium,  boron,  manganese,  gallium  and  rubidium).   The  focal 
point  of  the  problem  is  determined  to  be  the  recognition 
and  development  of  suitable  descriptors  to  be  used  in  the 
pattern  recognition  problem.   To  this  end,  outlines  for 
developing  numerical  estimates  of  such  descriptors  as 


IV 


linear  dynamic  range,  cost  of  analysis,  time  of  analysis 
and  selectivity  are  given.   Experimental  verification  of 
the  theoretical  definition  of  selectivity  is  given  by  the 
analysis  of  several  elements  in  a  simulated  sea  water  matrix 
via  atomic  absorption  and  atomic  emission  spectrometric 
techniques.   The  final  descriptor  proposed  is  a  theoretically 
calculated  limit  of  detection.   Although  perhaps  of  limited 
use  in  the  pattern  recognition  problem,  theoretically 
calculated  limits  of  detection  afford  the  analytical  chemist 
a  means  by  which  to  objectively  compare  ultimate  capabilities 
of  a  variety  of  spectroscopic  techniques.   Through  signal- 
to-noise  (S/N)  expressions,  limits  of  detection  are  calculated 
for  atomic  fluorescence  spectrometry,  atomic  absorption 
spectrometry,  atomic  emission  spectrometry,  optoacoustic 
spectroscopy  and  optogalvanic  spectroscopy.   After  identifying 
the  limiting  noise  sources  for  these  techniques,  the  S/N 

expressions  are  solved  for  the  limiting  observable  concen- 

3 
tration  (.in  atoms/cm  )  under  a  variety  of  experimental 

conditions.   Two-photon  processes  are  also  considered  (via 

virtual  and  real  levels)  for  atomic  fluorescence  spectrometry 

and  optogalvanic  spectroscopy.   Furthermore,  the  means  by 

which  one  can  convert  the  detection  limit  in  atoms/cm3  of 

hot  gases  to  ug/ml  Cppm)  of  analyte  in  solution  are  given. 

Comparisons  of  the  theoretical  limits  of  detection  to 

experimentally  derived  limits  of  detection  are  made  at  which 

time  good  agreement  between  the  two  is  observed. 


INTRODUCTION 

Analytical  chemistry  is  guided  on  its  course  to  ever 
increasing  sensitivity  and  selectivity  by  various  technolog- 
ical advances.   Accompanying  the  improvements  of  this  disci- 
pline as  a  science  are  the  rising  costs  of  instrumentation 
and  analysis.   No  longer  can  the  chemist  afford  himself  to 
take  too  many  "educated  guesses"  while  trying  to  solve  an 
analytical  problem.   It  would,  therefore,  be  of  great  eco- 
nomic importance  if  the  analytical  chemist  possessed  a 
systematic  and  objective  means  of  forecasting  the  utility 
of  a  particular  analytical  technique  to  solving  a  specific 
problem  at  hand.   In  general,  the  approach  taken  is  a 
purely  subjective  one  based  predominantly  on  the  personal 
experiences  of  the  analyst.   Although  perfectly  adequate 
at  times,  there  are  many  disadvantages  to  relying  solely 
upon  a  chemist's  intuition.   If  we  consider  the  rapidity 
with  which  new  techniques  and  methodologies  enter  the 
scientific  literature,  it  becomes  virtually  impossible  to 
keep  abreast  of  the  latest  advancements,  much  less  draw 
objective  conclusions  concerning  their  suitability.   The 
problems  associated  with  the  bringing  together  of  appropri- 
ate analytical  methodologies  with  analytical  problems 
can  be  classified  as  part  of  a  more  universal  concern, 


namely,  the  storage,  transf errence  and  utilization  of  large 
amounts  of  information.   To  this  end,  man  has  developed  the 
modern  day  computer.   However,  the  computer  is  rendered 
useless  unless  adequate  descriptors  of  the  information  to 
be  encoded  are  developed.   Furthermore,  algorithms  need  to 
be  developed  that  mimic  the  decision  making  process,  which 
will  in  turn  make  use  of  these  "subject ive-freeu  descriptors. 
The  ultimate  process  envisioned  is  one  in  which  the  analyti- 
cal chemist  can  define  the  problem  at  hand  by  the  utiliza- 
tion of  several  "key  words"  available  to  him,  and  subse- 
quently match  an  analytical  technique  to  suit  his  needs  via 
an  objective  mathematical  procedure. 

The  ensuing  work  has  been  performed  to  achieve  two 
purposes.   Primarily,  it  is  intended  to  show  that  the 
techniques  of  pattern  recognition  have  great  practical 
significance  in  solving  the  problems  of  analytical  chemistry. 
This  will  be  illustrated  by  a  novel  application  of  pattern 
recognition  to  the  discerning  of  the  geographical  origins 
of  frozen  concentrated  orange  juice  via  their  trace  element 
concentrations.   Once  the  viability  of  the  technique  has 
been  established,  a  discussion  of  its  application  to  bring- 
ing together  analytical  problems  and  techniques  will  follow. 
The  importance  of  establishing  impartial  and  sensitive 
descriptors  for  the  techniques  will  be  emphasized,  at 
which  time  the  proposed  descriptors  will  be  introduced; 
selectivity,  cost,  time  and  linear  dynamic  range  amongst 
others.   Emphasis  will  be  placed  on  making  all  proposed 
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figures  of  merit  arithmetically  quantifiable.   Attempts  will 
be  made  to  calculate  theoretical  limits  of  detection  for  a 
variety  of  atomic  spectroscopic  techniques  in  order  to 
establish  means  of  comparing  their  ultimate  capabilities. 
This  in  turn  may-  find  suitable  usage  in  the  pattern  recog- 
nition problem  as  well. 


CHAPTER  ONE 
PATTERN  RECOGNITION 

Introduction 
The  main  purpose  of  pattern  recognition  is  to  categorize 
a  sample  of  observed  data  as  a  member  of  a  class  to  which  it 
belongs.   Typical  examples  of  everyday  usages  of  pattern 
recognition  involve  the  evaluation  of  common  objects  such  as 
Toyotas,  German  Shepherds  and  Professors  as  part  of  a  larger 
class  of  automobiles,  canines  and  humans.   The  actual  mechanics 
of  pattern  recognition  were  not  investigated  until  the  1950' s 
with  the  advent  of  computers  and  the  birth  of  cognitive 
psychology  (1).   Since  then,  pattern  recognition  has  found 
many  and  varied  applications  which  are  dramatized  by  several 
excellent  reviews  which  have  appeared  in  the  literature 
(2-6).   The  more  notable  applications  of  pattern  recognition 
are  the  interpretation  of  printed  characters  (letters)  and 
script  (7),  speech  recognition,  fingerprint  analysis  (8), 
analysis  of  microphotographs  of  biological  material  including 
blood  cells  and  chromosomes  (9),  and  the  analysis  of  photo- 
graphs of  events  in  bubble,  spark,  or  cloud  chambers  (10). 
With  respect  to  chemical  applications,  pattern  recognition 
has  helped  organize  the  enormous  amount  of  complicated 


5 
spectra  generated  by  such  techniques   as  mass  spectrometry 

and  infrared  spectrometry.   It  has  helped  elucidate  structures 
of  compounds  when  combining  the  information  derived  from  both 
analysis. 

Before  discussing  the  present  example  of  the  application 
of  pattern  recognition,  it  would  be  instructive  to  acquaint 
the  reader  with  the  fundamental  mathematical  concepts  and 
terminology.   A  simple  physical  interpretation  will  serve  as 
a  convenient  starting  point  for  our  discussion.   Essentially, 
we  are  interested  in  describing  geometric  regions  in  space 
and  then  finding  a  means  to  discriminate  between  them.   The 
regions  in  space,  or  ''categories,"  are  described  by  "features^' 
or  "descriptors,"  of  the  measurement  process.   Each  descriptor 
is  associated  with  a  vector  whose  direction  defines  a  dimen- 
sion, or  axis  of  space,   and  whose  magnitude  defines  the 
distance  along  that  axis.   For  instance,  if  two  elements 
(A,B)   are  to  be  analyzed     quantitatively, in  a  particular 
sample  matrix,  we  will  be  dealing  with  a  2-dimensional  space 
whose  magnitude  along  any  axis  is  governed  by  the  measured 
concentration  for  that  element.   We  now  can  define  a  pattern 
vector  for  the  results  derived  from  one  analysis  of  the  two 
elements  of  interest  as 

X  =  (A,B).  (1) 

This  vector  now  represents  a  single  point  in  2-dimensional 
space.   If  we  now  make  the  assumption  that  samples  taken 
from  a  homogeneous  population  will  display  similar  concentra- 
tions for  the  2  test  elements,  then  their  associated  pattern 
vectors  will  define  points  in  a  similar  region  in  space. 
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We  now  have  a  basis  upon  which  to  form  a  category.   Further- 
more, with  a  second  population,  differing  from  the  first 
with  respect  to  the  concentrations  of  the  two  target  elements 
of  interest,  we  can,  by  analogy,  form  a  second  category. 
Having  formed  these  two  categories,  we  are  now  able  to  dis- 
tinguish between  the  two  sample  populations  by  merely  distin- 
guishing between  the  two  categories  formed  in  space.   By 
examining  Figure  1,  it  is  apparent  that  a  surface  can  be  drawn 
that  divides  the  two  categories.   This  now  transforms  the 
question  of  "to  what  category  does  a  point  belong?"  to 
"on  which  side  of  the  dividing  line  does  a  point  alleged  to 
be  in  one  of  the  categories  fall?1',.   This  dividing  line  is 
appropriately  termed  a  "decision  vector" 

D  =  (a,b)  (2) 

for  it  allows  us  to  make  our  final  categorization  decision. 
Upon  obtaining  the  dot  product  of  a  pattern  vector  with  a 
decision  vector,  the  pattern  vector  is  now  fixed  in  space 
with  respect  to  the  decision  vector.   In  the  ideal  case,  the 
dot  product  of  the  decision  vector  and  pattern  vector  will 
yield  one  of  two  results, 

X-D  =  {Aa  +  Bb}  =  S  (3) 

where  S>0  or  S<0.  Finally,  we  can  arbitrarily  define  the 
decision  vector  such  that  when  S>0,  then  X  will  belong  to 
category  I  and  when  S<0,  then  X  will  belong  to  category  II. 

It  should  be  pointed  out  at  this  time  that  rigorous 
mathematical  procedures  for  both  producing  pattern  vectors 
and  decision  vectors  have  been  formulated  (11-12).   The 
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Separation  of  Pattern  Vectors 
Derived  from  Two  Distinctly  Different  Populations 
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programs  developed  allow  patterns  to  be  formed  in  data  where 
upon  inspection  they  are  not  readily  apparent.   This  is 
extremely  helpful  in  making  otherwise  intractable  problems 
amenable  to  analysis  by  pattern  recognition.   However,  for 
the  less  complex  problems,  upon  careful  examination  of  the 
data,  an  appropriate  and  effective  decision  vector  can  be 
arrived  at  with  a  minimal  amount  of  effort.   In  the  end,  the 
efficiency  of  a  decision  vector  is  determined  by  repeated 
-attempts  at  making  the  discriminating  decision  and  recording 
the  success  rate. 

Discrimination  of  Geographical  Origins  of  Samples  of  Frozen 
Concentrated  Orange  Juice 

Having  covered  the  fundamentals  of  pattern  recognition, 
we  are  ready  to  show  its  application  to  the  present  work. 
The  objective  of  this  study  was  to  compare  a  number  of 
samples   of  Florida  and  Brazilian  orange,  jul^e  concentrates 
as  to  their  inorganic  elemental  content  for  the  purpose  of 
fingerprinting  their  geographical  origins  (13).   A  detailed 
sample  preparation  scheme  has  been  formulated  by  McHard 
et  al.  (.13)  which  essentially  entails  dry  ashing  20  g  samples 
of  the  frozen  concentrated  orange  juice  (FCOJ)  for  approxi- 
mately 16  hours.   Dissolution  of  the  ash  into  50  ml  of 
0.1  M  HNO3  is  then  followed  by  analysis  by  atomic  emission 
spectroscopy.   The  instrument  used  was  a  Spectrametrics ,  Inc., 
Spectraspan  II  equipped  with  the  Spectrajet  III  Plasma  source. 
Details  of  the  analysis  can  be  found  in  McHard  et  al .  (14), 
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while  a  compilation  of  their  results  can  be  found  in  Table  I. 
From  the  results, . McHard  et  al.  were  able  to  show  that  signif- 
icant differences  in  concentrations  existed  between  Floridian 
FCOJ's  and  Brazilian  FCOJ's  with  respect  to  barium,  boron, 
gallium,  manganese  and  rubidium  concentrations.   It  should 
be  mentioned  at  this  point  that  although  concentrations  in 
the  FCOJ  samples  should  ideally  be  analyzed  with  respect  to 
a  standard  reference  of  known  matrix  composition,  none  was 
readily  available  for  the  analysis.   Therefore,  concentrations 
are  reported  as  ratios  to  an  internal  reference  element,  zinc, 
which  was  shown  to  be  of  fairly  constant  concentration  in 
both  the  Floridian  and  Brazilian  populations.   Furthermore, 
the  use  of  a  particular  brand  of  Florida  juice  as  a  reference 
bypassed  the  need  to  account  for  solids  content  differences 
between  the  two  populations.   To  lend  statistical  credence 
to  the  conclusions  that  fingerprinting  elements  had  been 
established,  an  intensive  statistical  analysis  of  the  data 
was  performed.   By  use  of  Tschebychef f ' s  inequality,  numerical 
probabilities  can  be  assigned  to  the  chances  of  an  analyzed 
sample  being  of  Floridian  origin.   This  entailed  calculating 
means  and  standard  deviations  for  the  values  of  each  of  the 
elements  found  in  the  samples  analyzed  within  the  Florida 
population.   This  now  became  our  reference  population  with 
which  subsequent  unknown  analysis  could  be  compared.   We  then 
assign  a  deviation  'a'  for  a  test  sample,  which  is  a  measure 
of  the  number  of  standard  deviations  away  from  the  mean  of 
the  reference  population  that  this  value  lies.   From  this  a 
probability  (P)  of  the  test  sample  being  part  of  the  Floridian 
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population  can  be  calculated  using  Tschebychef f * s  inequality, 

P( |x  -  x| )  >  as  «  -   .  (4) 

This  method  is  quite  successful  at  achieving  its  intended 
goal   of  discriminating  between  the  two  populations.   However, 
regardless  of  its  intuitive  appeal  due  to  its  simplicity, 
the  calculations  needed  to  arrive  at  a  decision  are  tedious 
and  time  consuming  at  the  very  best.   What  is  proposed  is 
the  utilization  of  a  somewhat  less  common  mode  of  analysis 
of  the  data  which  also  has  great  discriminatory  power,  yet 
entails  a  much  less  laborious  manipulation  of  the  data.   The 
advantages  of  the  application  of  pattern  recognition  tech- 
niques to  this  problem  over  those  of  the  traditional  statis- 
tical methods  will  be  discussed  in  greater  detail  after 
the  problem  has  been  worked  out . 

The  same  five  elements  chosen  by  McHard  et  al.,  i.e. 
barium,  boron,  gallium,  manganese  and  rubidium,  are  the  ones 
chosen  for  this  analysis.   This  was  done  to  allow  a  direct 
comparison  between  the  two  methods.   In  actuality,  further 
optimization  of  the  pattern  recognition  technique  can  be 
realized  if  the  remaining  21  elements  were  objectively 
tested  (.11)  to  determine  their  suitability  as  discriminating 
descriptors.   However,  upon  inspection  of  the  data  base  of 
these  5  elements,  it  is  seen  that  they  exhibit  no,  or  minimal, 
overlap  in  their  ranges  which  does  indeed  allow  for  precise 
category  formation.   We  were  also  fortunate  enough  to  have 
the  absolute  values  of  the  ranges  of  concentrations  to  be 
approximately  equal  (see  Table  II).   If  this  were  not  so,  we 
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Table  II 


Summary  of  Concentration  Ranges  for  the 
Five  Fingerprinting  Elements  in 
Brazilian  and  Floridian  Samples 
of  Frozen  Concentrated  Orange  Juice 


Fla. CPpm)a 
Ba     0.025-0.07    x=0.048 
B     0.95-1.20     x=1.08 
Ga     0.03-0.04     x=0.035 
Mn     0.25-0.315    x=0.28 
Rb     0.365-0.74    x=0.55 


Brazil (ppm) 
0.139-0.53   x=0.33 
0.675-1.61   x=1.14 
0.047-0.11   x=0.079 
Q . 30-0 . 60    x=0 . 45 
1.98-4.86    x=3.42 


a  -  all  values  reported  are  adjusted  for  variations  in  solid 
content  and  express  concentrations  to  be  found  in  single 
strength  orange  juice 
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would  perceive  a  disproportionate  weighting  factor  assigned 
to  that  element  whose  concentration  values  were  higher  than 
the  remaining  ones.   However,  had  the  situation  arisen  we 
would  have  been  able  to  linearize  the  data,  thereby   pre- 
cluding this  possibility  (.12). 

Having  chosen  these  5  elements  as  our  descriptors,  we 
have  established  the  fact  that  we  will  be  forming  categories 
in  5-dimensional  space.   Obviously,  this  problem  cannot  be 
handled  graphically  nor  imagined  by  the  human  mind.   This 
illustrates  one  of  the  main  advantages  of  pattern  recognition: 
the  ability  to  handle  multidimensional  problems  as  easily  as 
it  handles  those  of  lesser  dimensions.   As  previously  men- 
tioned, the  magnitude  along  any  one  particular  axis  is  gov- 
erned by  the  measured  concentration  for  that  particular 
element.   We  now  can  define  a  pattern  vector  for  one  sample 
set  (one  particular  FCOJ  sample) 

X  -  (IBal,  [Bl,  [Ga3,  [Mnl,  [Rbl)  (.5) 

This  vector  now  represents  a  single  point  in  5-dimensional 
space  that  completely  characterizes  that  one  particular 
FCOJ  sample.   If  we  now  make  the  assumption  that  similar 
FCOJ  samples  can  be  described  by  pattern  vectors  that  define 
points  in  a  similar  region  in  space,  we  now  have  a  basis 
upon  which  to  form  categories.   The  problem  now  reduces  to 
one  in  which  we  are  interested  in  discriminating  between 
these  categories.   In  practice,  we  are  looking  for  a  means 
whereby  we  can  classify  a  sample  into  one  of  two  categories, 
i.e.,  Brazilian  FCOJ  or  Floridian  FCOJ.   The  means  by  which 
this  is  accomplished  is  by  the  introduction  of  the 
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aforementioned  decision  vector .  After  .careful  inspection  of  the 
data  in  Table  II,  the  following  decision  vector  was  arrived 
at  j 

D  =  CI,  -1,  -1,  -1,  1)  (6) 

The  decision  vector  was  formulated  such  that  if  the  dot 
product  of  it  and  a  sample  vector  (,S)  is  less  than  zero, 
then  X  is  a  Floridian  juice,  and  if  S  is  greater  than  zero, 
then  X  is  a  Brazilian  juice.   The  fact  that  the  dot  product 
would  result  in  either"  a  positive  or  negative  integer  is 
just  a  matter  of  convenience.   Any  arbitrarily  chosen  de- 
cision vector  would  be  sufficient  in  defining  the  area  in 
space  that  contains  the  categories  of  interest.   Their 
boundaries,  however,  would  not  be  as  conveniently  defined 
as  in  the  present  example.   Nevertheless,  once  the  boundaries 
are  evaluated  they  can  be  utilized  as  equally  effective 
vehicles  of  the  decision  making  process. 

In  general,  trace  elemental  concentrations  in  juices 
of  Floridian  origin  were  lower  than  those  of  non-domestic 
origin.   Therefore,  three  cases  were  examined  to  detect 
the  efficiency  of  the  decision  vector. 

Case  I .   Here,  we  wish  to  distinguish  between  a 
Florida  juice  with  the  lowest  possible  values  and  a 
Brazil  juice  with  the  highest  possible  values. 

Case  II.   Here,  we  wish  to  distinguish  between  a 
Florida  juice  with  average  values  and  a  Brazil  juice  with 
average  values. 
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Case  III.   Here,  we  wish  to  distinguish  between  a 
Florida  juice  with  the  highest  possible  values  and  a  Brazil 
juice  with  the  lowest  values.   Case  III  is  the  most  important 
case  because  it  represents  the  case  where  maximal  overlap 
between  the  two  origins  may  exist.   The  concentrations  for 
any  of  the  five  target  elements  are  assumed  to  be  uncorre- 
cted to  one  another.   If  this  assumption  is  true,  and  ex- 
perimental evidence  indicates  this,  then  the  probability 
of  finding  a  sample  that  has  the  concentrations  of  all 
five  of  the  target  elements  stacked  to  one  side  of  the  mean 
is  nominal.   Therefore,  case  III  actually  represents  the 
worst  possible  case  we  can  ever  expect  to  encounter  and 
does  not  portray  the  average  real  sample.   It  is,  however, 
the  ultimate  limiting  case  and  will  serve  as  a  good  indi- 
cator of  pattern  recognition's  ability  to  discriminate 
between  the  two  geographical  origins.   Referring  to  Table 
II,  cases  I-III  will  be  evaluated: 

Case  I.   a)  X1   -  (.0.025,  0.95,  0.03,  0.25,  0.37) 

(Floridian) 

where  the  above  vector  represents  the  following  ordering  of 

elements: 

(Ba,  B,  Ga,  Mn ,  Rb) 

and  the  decision  vector  was  given  by  equation  (6)  as 

D  =  (1,  -1,  -1,  -1,  1). 

We  must  evaluate  the  dot  product  of  the  two  vectors,  D-X^, 

D-X-l  =  (.0.025  -  0.95  -  0.03  -  0.25  +  0.37)  =  S 

S  =  -0.84 
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S<0;  therefore,  X-j_  is  a  Floridian  juice. 

b)  X2  =  CO. 53,  1.61,  0.11,  0.60,  4.86) 

(.Brazilian) 

D-X2   =  3.07 

S>0;  therefore ,  X2  is  a  Brazilian  juice. 


Case  II.   a)  X±  =    (0.048,  1.08,  0.035,  0.28,  0.55) 

(Floridian) 

D-Xi  =  -0.80 

S<0;  therefore,  X-^  is  a  Floridian  juice. 

b)  X2  =  (0.33,  1.14,  0.079,  0.45,  3.42) 

(Brazilian) 

D-X2  =  2.49 

S>0;  therefore,  X2  is  a  Brazilian  juice. 


Case  III.   a)  X±   =  (.0.07,  1.20,  0.04,  0.315,  0.740) 

(Floridian) 

D-Xi  ■  -0.745 

S<0; therefore,  X1  is  a  Floridian  juice. 

b)  X2  =  (0.139,  0.675,  0.047,  0.30,  1.98) 


D-X2  =  1.10 


(Brazilian) 


S>0; therefore,  X2  is  a  Brazilian  juice. 

The  results  of  the  analysis  of  the  five  target  elements 
were  in  actuality  ratioed  to  zinc.   To  make  the  procedure 
functional  for  this  set  of  data,  the  decision  vector  had  to 
be  modified  to  reflect  the  ratios  as  tabulated  in  Table  I. 
The  new  decision  vector  is 

Dn  =  CI,  -2,  -2,  -2,  1). 
We  will  now  repeat  case  III  to  substantiate  the  efficiency 
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of  the  new  decision  vector. 


Case  III.   a)  X±   =  (2.9661,  0.8702,  1.2812,  1.0586, 

1.8473) 

(JFloridian) 

Dn-X1  =  -1.603 

S<0;  therefore,  X-^   is  a  Floridian  juice. 

b)  X2  =  (.5.150,  0.9448,  1.3794,  1.3816,  6.0511) 

(Brazilian). 

Dn-X2  =  3.789 

S>0; therefore,  X2  is  a  Brazilian  juice. 


Conclusion 
It  has  been  shown  that  for  the  worst  possible  case  in 
which  all  elements  of  interest  in  a  Florida  FCOJ  have  their 
highest  possible  values,  and  those  of  the  Brazilian  FCOJ 
have  their  lowest ,  the  proposed  technique  offers  a  clear 
and  substantial  separation  of  the  two  groups.   With  the 
sample  of  the  population  that  has  been  obtained  to  date, 
the  decision  vector  is  100%  successful.   Furthermore, 
after  having  recognized  that  patterns  were  being  formed  by 
the  data,  the  application  of  pattern  recognition  techniques 
to  the  problem  proved  trivial.   This  is  a  point  which  the 
author  wishes  to  emphasize.   Although  it  is  recognized 
that  pattern  recognition  can  handle  problems  that  otherwise 
are  too  complex,  it  also  is  a  viable  technique  to  use  when 
conventional  statistical  and  probabilistic  methods  can 
ordinarily  be  used.   The  distinct  advantage  comes  from  the 
time  saving  aspect  of  the  minimized  number  of  calculation 
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steps  needed  in  order  to  arrive  at  a  decision.   Once  the  data 
base  is  established  and  a  proper  decision  vector  is  formed, 
the  actual  testing  of  the  data  becomes  very  rapid.   It  is 
the  author's  opinion  that  this  alone  makes  the  pattern 
recognition  technique  attractive.   There  is  one  other  ad- 
vantage that  may  not  be  so  apparent  and  should  be  pointed 
out.   This  lies  in  the  fact  that  once  a  person  is  exposed 
to  the  techniques  of  pattern  recognition,  he  begins  to 
cognizantly  think  in  terms  of  groups  or  patterns.   This 
allows  interrelationships  between  items  to  be  scrutinized 
from  a  natural  perspective,  i.e.  in  terms  of  patterns. 
Once  exposed  to  the  techniques  of  pattern  recognition,  and 
after  familiarizing  oneself  with  the  basic  mathematical 
tools  of  the  technique,  more  and  more  problems  will  become 
amenable  to  being  tackled  by  pattern  recognition.   Pattern 
recognition  is  a  powerful  mathematical  technique  that 
allows  strong  conclusions  to  be  drawn  about  sets  of  data 
with  great  confidence.   It  may  not  replace  more  traditional 
statistical  methods,  but  can  definitely  serve  as  an  attrac- 
tive alternative. 


CHAPTER  TWO 
METHOD  COMPARISON 

Introduction 

As  was  stated  in  the  introduction,  one  of  the  major  goals 
of  this  work  is  to  show  the  applicability  of  pattern  recog- 
nition  techniques  to  the  problem  of  choosing  between  a  variety 
of  analytical   tecniques  to  solve  a  specific  analytical 
problem.   Many  authors  have  addressed  themselves  to  the  problem 
of  objective  method  comparison  and  a  variety  of  approaches 
have  been  developed. 

One  of  the  simplest  and  most  straightforward  approaches 

to  method  comparison  was  first  presented  by  McFarren  et  al. 

(15).   Their  approach  was  based  on  the  assumption  that  the 

only  criterion  of  interest  to  an  analyst  about  to  chose 

between  several  analytical  techniques  to  solve  a  particular 

problem  is  its  error  content.   A  proposed  "total  error" 

of  an  analytical  technique  is  calculated  and  criterion  by 

which  decisions  will  be  made  based  on  the  total  error 

calculations  is  established.   The  formula  for  calculating 

the  total  error  of  an  analytical  technique  is 

|E  |  +  2(s) 

— - — = X  100  =  Total  error  (1) 

Xt 
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where ■ 

E   =  absolute  value  of  mean  error 
m 

s   =  standard  deviation  of  the  measurements 

X  =  true  value  of  the  analyzed  component . 
The  mean  error  is  defined  as  the  difference  between  the  average 
of  a  series  of  test  results  and  the  true  result.   The  standard 
deviation  retains  its  classical  definition  which  can  be  found 
in  any  elementary  statistics  book. 

Upon  calculation  of  the  total  error,  which  is  actually 
a  measure  of  the  precision  plus  accuracy  of  the  technique,  it 
remains  to  be  compared  to  the  accepted  evaluation  criterion 
and  a  decision  is  then  made.   The  proposed  criterion  was  that 
any  technique  that  exhibited  a  total  error  of  less  than  25% 
was  judged  to  be  excellent.   Those  techniques  with  a  total 
error  greater  than  25%  but  less  than  50%  were  judged  to  be 
acceptable  while  those  having  total  errors  greater  than  50% 
were  judged  unacceptable.   In  order  to  make  this  approach 
viable,  certain  variables  have  to  be  strictly  defined.   For 
instance,  methodologies  for  sample  analysis  have  to  be 
optimized  and  explicitly  recorded.   This  is  true  under  any 
and  all  circumstances  if  meaningful  comparisons  are  to  be 
made.   Instruments  must  be  carefully  calibrated,  reference 
standards  must  be  established  and  proper  blank  correction 
must  be  performed.   There  are  several  major  drawbacks  to  this 
approach  which  should  be  pointed  out.   Primarily,  the  total 
error  of  analysis  may  not  be  the  sole  determining  factor  of 
the  applicability  of  an  analytical  technique.   The  cost  of 
the  analysis  is  always  a  major  consideration  when  faced  with 
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various  budgetary  constraints.   Can  the  same  information,  or 

at  least  information  that  will  still  suit  the  needs  of  the 
analyst,  be  obtained  at  a  lower  cost?   Another  factor  that 
must  be  considered  is  the  speed  of  the  analysis.   It  would  be 
foolish  to  choose  a  particular  analytical  method  that  required 
time  consuming  preparative  steps  when  a  direct  analysis  can 
be  done  (assuming  all  other  variables  remain  constant).   What 
all  this  is  leading  to  is.  that, when  faced  with  making  the 
type  of  decision  we  are  discussing  here,  all  pertinent  figures 
of  merit  for  the  proposed  analytical  techniques  must  be  looked 
at.   It  is  insufficient  to  just  consider  the  precision  and 
accuracy  obtainable  by  any  particular  analytical  method.   One 
of  the  first  scientists  to  identify  the  need  to  establish  good 
qualitative  figures  of  merit  which  are  quantifiable  was 
Heinrich  Kaiser  (16).   It  was  his  goal  to  reduce  the  task  of 
choosing  between  analytical  methods  to  a  purely  mathematical 
process.   Using  the  ideas  and  terminology  of  pattern  recognition, 
one  would  have  to  represent  the  different  possible  analytical 
problems  as  points  in  a  multi-dimensional  space.   With  this 
designation,  related  problems  would  inevitably  become  neigh- 
boring points  in  this  space.   Furthermore,  a  corresponding 
order  must  be  given  to  the  various  analytical  procedures, 
which  can  be  represented  as  points  in  the  space  of  procedures. 
The  problem  now  becomes  reduced  to  one  of  establishing  a 
correlation  between  points,  or  regions ,  in  the  two  spaces. 
Or,  if  both  procedures  and  problems  are  defined  by  the  same 
descriptor,  and  we  are  interested  in  determining  what  procedure 
to  use  for  a  particular  problem,  the  problem  is  reduced  to  one 
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of  determining  the  appropriate  category  to  which  the  problem 
descriptors  belong.   This  is  analagous  to  the  problem  of 
deciding  upon  the  geographical  origins  of  a  particular 
frozen  concentrated  orange  juice  sample  which  was  discussed 
previously. 

As  with  any  pattern  recognition  problem,  the  effect- 
iveness of  its  application  is  solely  determined  by  the  chosen 
descriptors.   They  must  be  good  indicators  of  the  techniques 
and  should,  ideally,  be  independent  of  one  another.   A  more 
fundamental  necessary  condition  for  any  effective  descriptor 
is  its  ability  to  be  assigned  discrete  numerical  values.   It 
should  be  fairly  apparent  that  when  trying  to  associate  a 
particular  analytical  method  with  a  specific  problem,  it  is 
important  to  know  the  limit  of  detection  of  the  techinique, 
its  linear  dynamic  range,  its  selectivity  and  specificity, 
the  approximate  speed  and  cost  involved,  the  sample  size 
needed  and  its  allowable  physical  and  chemical  form,  the 
accuracy  and  precision  that  is  obtainable,  as  well  as  whether 
or  not  the  instrument  is  available  to  the  analyst.   All  of 
these  figures  of  merit  seem  intuitively  appropriate  to  the 
method  comparison  problem;  however,  it  is  not  readily  apparent 
how  to  quantify  some  of  them. 

To  date,  only  one  attempt  has  been  made  to  actually 
apply  the  techniques  of  pattern  recognition  to  the  problem  of 
deciding  amongst  a  variety  of  analytical  techniques  to  solve 
a  particular  analytical  problem  (17).   This  work  concerned 
itself  with  discriminating  between  the  applicability  of 
UV-Vis  Molecular  Absorption  Spectrophotometry  (MA)  and 
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Atomic  Absorption  (AA)  Spectrophotometry  to  specific  analytical 
problems.   The  author  defined  several  descriptors  and  was 
able  to  assign  quantitative  ranges  for  their  values.   The 
questions  Dr.  Vandeginste  asked  whose  answers  defined  the 
descriptors  used  are  as  follows: 

a)  Is  the  analysis  for  qualitative  or  quantitative  work? 

b)  What  is  the  accuracy? 

c)  What  is  the  cost  of  the  analysis? 

d)  What  is  the  analysis  time? 

e)  Is  the  analysis  destructive  or  nondestructive? 

f)  What  is  the  sample  quantity  needed  for  analysis? 

g)  What  is  the  concentration  to  be  determined? 
h)  What  is  the  compound? 

i)  What  type  of  matrix  is  it  in? 
After  surveying  the  literature  for  problems  that  have  already 
been  solved  by  either  MA  or  AA,  a  total  of  237  articles  were 
chosen.   Of  these  237,  157  were  chosen  to  produce  the 
methodology  space.   The  remaining  80  papers  were  to  be  used 
as  the  problem  vectors  and  were  to  be  categorized  with 
respect  to  the  two  categories  formed  in  methodology  space. 
The  results  of  this  approach  were  surprisingly  good  consider- 
ing the  fact  that  the  chosen  figures  of  merit  are  not 
necessarily  the  optimal  ones  to  be  used.   Correct  classifica- 
tions were  made  on  the  order  of  85%  of  the  time.   The  major 
obstacle  to  Dr.  Vandeginste ' s  work,  as  well  as  to  that  of 
any  researcher  in  the  field,  is  the  lack  of  proper  documenta- 
tion found  in  the  literature.   For  instance,  although  the 
cost  of  an  analysis  would  seem  to  be  a  figure  of  merit  when 
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comparing  analytical  methodologies,  very  few  authors,  choose 
to  report  the  dollar  and  cents  figures  involved  with  their 
analysis.   Furthermore,  the  time  of  an  analysis  seems  as  if 
it  should  be  a  crucial  deciding  factor  when  chosing  amongst 
analytical  techniques,  yet,  it  too  is  seldom  reported  in  the 
literature.   The  only  categorical  exception  to  this  is  the 
reporting  of  elution  times  in  chromatographic  analyses. 
However,  even  in  chromatographic  methods,  elution  times  may 
not  be  the  sole  determinants  of  the  total  analysis  time. 
There  may  be  various  steps  in  sample  pretreatment  which  may 
or  may  not  be  the  rate  determining  steps  with  respect  to  the 
entire  analytical  procedure.   The  omission  of  these  and  other 
figures  of  merit  from  the  literature  is  not  a  mere  oversight 
of  the  respective  authors.   More  precisely,  it  appears  as  if 
the  need  to  report  these  figures  of  merit  has  not  been 
established.   This  perhaps  stems  from  the  lack  of  a  clear 
and  concise  format  for  reporting  them.  To  illustrate  some 
of  the  difficulties  involved  with  precisely  defining  a 
figure  of  merit  and  assigning  it  a  numerical  value,  let  us 
look  at  a  few  examples. 

Figures  of  Merit 
Time  of  Analysis 

If  a  particular  laboratory  worker  is  interested  in 
adopting  a  new  analytical  method  to  analyze  for  a  particular 
sample  and  is  particularly  interested  in  minimizing  the  time 
of  analysis,  assuming  all  other  figures  of  merit  remain 
equal,  how  should  he  go  about  doing  this?   To  begin  with, 
before  actually  performing  the  analysis  a  technician  will  have 
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to  be  trained  on  the  particular  instrument  involved.   The 
time  taken  to  accomplish  this  may  vary  drastically  depending 
on  the  experience  of  the  analyst  and  the  complexity  of  the 
analysis.   This  is  an  important  consideration  with  respect 
to  the  total  analysis  time  but  will  not  be  quantified  due  to 
its  great  variability.   Assuming  that  adequate  training  has 
been  provided  we  now  can  take  a  closer  look  at  the  actual 
analysis  time  involved.   It  should  be  kept  in  mind  that  we 
are  hoping  to  achieve  a  maximum  amount  of  uniformity  and  the 
ability  to  define  precisely  the  total  analysis  time  such  that 
interlaboratory  comparisons  could  be  monitored  objectively. 
In  the  following  discusssion,  we  will  be  considering  the 
analysis  of  a  sample  by  instrumental  means  only;  however, 
a  more  general  approach  can  be  realized  to  encompass  any  mode 
of  analysis. 

Standard  and  reagent  preparation.   If  we  designate  the 
total  amount  of  time  needed  to  prepare  stock  solutions  of 
standards  and  necessary  reagents  as  (t),  and  the  number  of 
analysis  that  these  stock  solutions  will  last  for  as  (n),  then 
we  can  define  the  total  time  needed  to  prepare  all  necessary 
solutions  per  analysis  as 

t   =  t/n      (minutes/analysis) 

Pretreatment .   Any  time  necessary  to  pretreat  the  sample 
and  make  it  amenable  for  analysis  must  be  considered.   We 
will  designate  this  time  interval  as  (t  )  and  will  include 
the  time  needed  for  any  physical  or  chemical  treatments  of 
the  samples  as  needed. 
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Analysis.   With  respect  to  the  instrument  itself,  we 
assume  that  suitable  warmup  to  attain  equilibrium  can  be 
achieved  without  operator  attention.   Therefore,  the  time 
required  for  analysis  is  defined  as  the  time  interval  between 
the  introduction  of  the  first  sample  or  standard  into  the 
instrument  and  the  time  it  takes  to  receive  adequate  infor- 
mation from  the  instrument  to  allow  calculation  of  the  desired 
quantity.   This  time  span  can  be  designated  by  (t  ). 

Calculation.   Finally,  once  adequate  information  in  the 

form  of  voltage  measurements,    digital  readouts, et c. ,  is 

obtained,  the  time  needed  to  convert  them  into  their  final 

usable  form  is  designated  (t  ).   If  the  instrument  is 

c 

interfaced  with  a  computer  this  time  span  is  nominal. 
However,  in  the  case  where  it  is  not,  this  time  interval  may 
become  significant. 

To  arrive  at  the  total  analysis  time  the  four  components 
are  added  algebraically, 

t,  ,  ,  =t   +  t   +t   +t 
total    r    p    a    c 

We  now  have  a  somewhat  unperturbing  method  to  calculate  a 
total  analysis  time  which  then  can  be  introduced  into  the 
pattern  recognition  problem.   The  effort  involved  to 
calculate  the  total  time  of  analysis  is  minimal  on  the  part 
of  the  analyst,  yet,  the  results  can  prove  very  beneficial 
to  method  comparison  studies.   Hopefully,  with  explicit 
directions  on  how  to  go  about  performing  these  calculations, 
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cooperating  researchers  would  make  the  effort  to  systeraize 
the  reporting  of  their  results  in  the  literature. 

Cost  of  Analysis 

Similarly,  a  systematic  approach  to  calculating  a  total 
cost  of  analysis  can  be  developed  (18,19).   The  importance 
of  this  figure  of  merit  can  easily  be  appreciated  and  a 
straightfoward  method  for  calculating  a  rough  estimate  would 
be  welcomed.   The  major  considerations  are  fourfold: 

1  -  cost  of  labor,  C, 

2  -  cost  of  instrument,  C. 

'   l 

3  -  cost  of  instrument  maintenance,  C 

'   m 

4  -  cost  of  supplies,  C  . 

To  estimate  the  cost  of  labor  required  to  perform  an 
analysis,  one  must  simply  consider  the  total  analysis  time 
(ttotal)  and  multiply  this  value  by  the  technician's  hourly 
wages.   The  initial  cost  of  the  instrument  should  be  available 
and  after  including  the  allowable  depreciation  due  to  wear 
and  tear,  the  total  yearly  cost  could  be  calculated  if  a 
finite  expected  lifetime  is  estimated.   If  records  are 
available  concerning  the  number  of  analyses  performed  per 
year  then  it  is  trivial  to  calculate  the  instrumental 
capital  cost  per  analysis.   With  respect  to  the  maintenance 
of  the  instrument,  careful  records  should  be  kept  detailing 
the  number  of  service  calls  needed  to  keep  the  instrument  in 
operating  condition.   If  repairs  are  performed  in  house  a 
logging  of  the  man  hours  spent  in  repairs  is  necessary. 
Calculating  the  cost  of  maintenance  per  analysis  requires 


29 

knowledge  of  the  average  cost  spent  on  a  yearly  basis  and 

dividing  by  the  number  of  analyses  performed  during  the  same 

time  period.   This  should  include  estimates  of  replacement 

parts  as  well.   Estimating  the  amount  of  money  spent  on 

supplies  may  be  a  little  more  difficult  than  the  preceding 

calculations.   The  majority  of  items  used  in  an  analysis  such 

as  glassware,  chemicals,  etc. ,  need  to  be  considered  but  could 

be  included  in  general  overhead  incurred  while  running  a 

laboratory.   The  only  upkeep  expenses  that  need  be  considered 

are  the  cost  of  specialty  gases  used,  power  consumption,  etc. 

For  the  most,  part  the  price  for  this  aspect  of  the  analysis 

with  respect  to  the  entire  analytical  procedure  is  negligible 

and  will  not  be  considered.   Totaling  the  individual  cost 

factors  mentioned  should  give  a  reasonable  estimate  of  the 

total  cost  of  an  analysis, 

C+n+al    =  C,  +  C.  +  C   +  C  . 
total    1    i    m    s 

In  most  cases  if  a  specific  task  is  to  be  accomplished  the 
cost  per  analysis  may  be  trivial  with  respect  to  the  value 
of  the  information  obtained  (20).   However,  to  treat  the 
problem  objectively,  thereby  enabling  an  unbiased  decision 
to  be  arrived  at  concerning  methodology  choices,  the  cost  per 
analysis  should  be  included  as  a  figure  of  merit.  Depending 
on  the  relative  emphasis  placed  on  any  one  figure  of  merit 
being  utilized  in  a  pattern  recognition  problem,  a  relative 
weighting  factor  can  be  assigned  (21),  thereby  stressing  the 
most  important  criteria  to  be  used  and  slighting  those  of 
lesser  importance. 
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Linear  Dynamic  Range 

Continuing  on  our  quest  to  define  appropriate  figures 
of  merit  for  analytical  techniques,  we  now  will  consider  the 
linear  dynamic  range  associated  with  the  response  of  the 
analytical  method.   Although  it  is  obvious  that  knowledge 
of  such  a  "performance  characteristic''  (22)  would  help  define 
the  extent  of  applicability  of  an  analytical  method,  here 
too  we  are  faced  with  a  lack  of  precise  definitions  in  the 
literature  (23). 

With  respect  to  the  upper  limit  of  the  concentration 

range,  regardless  of  how  it  is  defined  it  may  usually  be 

increased  by   (1)  diluting  the  sample,  (2)  using  a  less 

sensitive  response  parameter  (e.g.,  a  less  intense  spectral 

line,  etc.  ),  or   (3)  altering  the  instrumental  configuration 

(24).   If  any  of  the  above,  are  done  to  alter  the  upper  limit 

of  the  concentration  range  than  a  new  analytical  method  is 

produced  by  the  strictest  definition  of  "analytical  method," 

An  analytical  method  is  to  be  regarded  as  the  set  of 
written  instructions  completely  defining  the  procedure 
to  be  adopted  by  the  analyst  in  order  to  obtain  the 
required  analytical  result   (25): 

Therefore,  if  these  alterations  are  made  they  should  be 

included  in  the  description  of  the  analytical  method. 

The  upper  limit  of  concentration  can  be  defined  in  a 

variety  of  ways.   For  instance,  it  can  be  defined  as  that 

concentration  which  represents  the  point  in  the  calibration 

curve  immediately  preceding  a  deviation  from  the  theoretical 

relationship  between  concentration  and  intensity.   This, 

however,  would  be  inadequate  in  that  linear  calibration 

curves  are  not  mandatory  for  performing  quantitative  analysis. 
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Another  possible  definition  would  be  one  that  defined  the 
upper  limit  of  concentration  to  be  the  highest  concentration 
used  to  prepare  the  calibration  curve.   This  is  somewhat  of 
a  better  definition  in  that  it  specifies  the  limits  to  which 
the  analytical  technique  has  been  tested,  although  there  is 
a  large  amount  of  arbitrariness  associated  with  picking  the 
concentrations  used.   A  more  specific  and  useful  definition 
would  be  one  that  chooses  the  upper  limit  of  concentration 
as  the  greatest-  concentration  for  which  quantitative  tests 
of  the  precision  and  bias  of  analytical  results  have  been 
made.   This  definition  actually  expresses  the  greatest 
concentration  that  has  undergone  experimental  scrutiny  and 
whose  results  for  accuracy  and  precision  are  available  for 
inspection.   However,  it  must  be  stressed  that  although  this 
definition  appears  to  be  the  best  one  for  the  upper  limit  of 
concentration  its  value  can  easily  be  experimentally  varied. 

Conversely,  a  more  invariable  value  for  concentration  is 
obtained  when  considering  the  lower  limit  of  the  concentration 
range.   In  this  case,  the  limit  of  detection  of  the  analytical 
technique  is  usually  firm  and  can  be  improved  only  after 
considerable  instrumental  optimization.   Although  there  seems 
to  be  general  agreement  as  to  how  to  define  the  detection 
limit  of  an  analytical  method  today,  it  appears  as  if  there 
has  not  always  been  a  clear  statistical  definition  available. 
Excellent  review  articles  concerned  with  the  historical- 
development  of  the  concepts  behind  defining  a  limit  of 
detection  have  been  presented  by  Kaiser  (26,  27)  and 
Currie  (28).   The  accepted  definition  of  the  limit  of 
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detection  is  (29)  that  concentration  that  would  give  a 
signal  that  is  equivalent  to  three  times  the  standard  deviation 
of  the  blank  measure.   Algebraically  this  can  be  expressed 
as  follows: 


C,  ■   =  3(or   )  (S) 

lim    v  bl 


where 


C,  .   =  concentration  at  detection  limit 
lim 

abl  =  s"tandard  deviation  of  the  blank 
S  =  slope  of  the  calibration  curve. 
Having  defined  the  upper  limit  and  lower  limits  of  the 
concentration  ranges  that  are  produced  by  analytical  methods, 
we  now  have  a  method  for  quantitating  their  values.   These 
values  can  then  be  easily  incorporated  into  the  pattern 
recognition  problem. 

Selectivity 

Another  very  important  figure  of  merit  for  any  analytical 

procedure  is  its  degree  of  selectivity  for  the  analyte  of 

interest  over  any  interfering  constituent  of  the  sample  or 

blank  matrix.   In  1933,  G.  E.  F.  Lundell  published  an  article 

in  which  the  goals  of  the  analytical  chemist's  work  in  terms 

of  developing  selective  techniques  were  set  forth, 

The  determinator ' s  salvation  lies  in  the  development 
of  selective  methods  of  analysis,  and  his  final  resting 
place  will  be  a  heaven  in  which  he  has  a  shelf  containing 
92  reagents,  one  for  each  element,  where  No.  13  is  the 
infallible  specific  for  Al ,  No.  26  the  sure  shot  for  Fe, 
No.  39  the  unfailing  relief  for  Y,  and  so  on  to  U   (30). 

Not  only  has  Lundell' s  vision  been  made  more  inaccessible 

by  the  introduction  of  at  least  13  more  elements  in  the  past 
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50  years,  but  his  hopes  have  never  come  close  to  being  realized. 

Nevertheless,  his  implication  is  that  no  analytical  method 

should  be  judged  acceptable  unless  a  measure  of  its  selectivity 

has  been  performed.   Lundell  goes  on  to  say, 

There  is  no  dearth  of  methods  that  are  entirely 
satisfactory  for  the  determination  of  elements  when 
they  occur  alone.   The  rub  comes  in  because  elements 
never  occur  alone,  for  nature  and  man  frown  upon 
celibacy.   Methods  of  determination  must  therefore 
be  judged  by  their  selectivity.   It  is  in  this  respect 
that  most  methods  are  weak  and  improvements  must 
come   (30). 

It  is  true  that  experimental  improvements  must  be  ad- 
dressed to  enhancing  selectivity  as  well  as  limit  of  detect- 
ability,-  however,  by  what  criterion  does  one  judge  selectivity? 
In  the  colloquial  sense, the  meaning  of  selectivity  is  obvious 
but  how  does  one  go  about  defining  it  quantitatively?   In  the 
areas  of  ion  selective  electrodes  and  ion  chromatography, 
efforts  have  been  made  to  define  selectivity  coefficients. 
However,  no  general  approach  to  experimentally  defining  a 
selectivity  for  an  analytical  procedure  has  been  developed. 
The  theoretical  groundwork  defining  selectivities  has  been 
formulated  by  Heinrich  Kaiser  (16);  however,  prior  to  this 
work, no  experimental  verification  of  its  validity  has  been 
reported.   If  a  numerical  analysis  of  the  selectivity  of 
an  analytical  procedure  can  be  developed  an  invaluable  and 
sensitive  figure  of  merit  will  result  that  can  effectively 
be  introduced  into  the  pattern  recognition  problem.   The 
following  work  has  been  undertaken  with  just  such  a  purpose 
in  mind. 

An  analytical  procedure  is  said  to  be  fully  ''selective'-' 
when  it  enables  the  different  components  in  the  sample  to  be 
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determined  independently  of  each  other.   Therefore,  the 
analyst  using  such  a  procedure   can  select  which  of  the 
components  are  to  be  determined  and  which  are  not,  without 
any  risk  or  errors  due  to  disregarding  some  components  of 
minor  interest.   The  phrase  ''more  or  less  selective"  indicates 
that  selectivity  is  understood  to  be  a  metric  quantity  and 
Kaiser  endeavoured  to  define  it  mathematically.   This  would 
enable  numerical  comparisons  to  be  made  of  the  varying 
degrees  of  selectivity  of  analytical  techniques  in  general. 

Theory 

The  numerical  value  of  selectivity  (E)  of  a  given 
analytical  procedure  can  be  derived  from  the  system  of  cali- 
bration functions  which  represent  the  relationship  between 
concentration  (c)  and  measurable  quantities  (x)  in  a  multi- 
component  analysis.   A  general  derivation  has  been  presented 
by  Kaiser  (31)  but  for  the  special  case  of  trace  elemental 
analysis,  which  we  will  consider  here,  a  simpler  method  leads 
to  the  same  results. 

Consider  a  multicomponent  system  comprised  of  three 
components  to  be  analyzed  spectroscopically  for  their  trace 
elemental  content.   If  the  elements  exhibit  no  spectral 
overlap  within  the  spectral  bandpass  of  the  monochromator  and 
if  no  interelement  matrix  effects  are  present,  then  the 
analysis  can  be  performed  in  a  straightf oward  manner.   The 
analytical  signals  measured  are  functions  of  one  component 
at  any  one  particular  wavelength  and  no  deconvoluting  of 
interferences  has  to  be  performed.   However,  if  spectral 
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overlaps  do  exist,  then  the  measured  intensity  is  a  function 
of  all  three  components  at  any  one  wavelength.   The  contri- 
bution to  the  total  signal  by  the  analyte  has_to.be.  ... 
determined  before  any  quantitative  information  about  the 
analyte  can  be  obtained.   Experimentally,  this  is  done  by 
measuring  the  intensity  of  the  three  component  system  at  the 
three  analytical  wavelengths f thereby  developing  a  series  of 
three  simultaneous  equations  called  the  calibration  functions. 
These  are  represented  as  follows, 

xl  =  H(C1>    C2>    C3} 

x2  =  g2(c1,  c2,  c3)  (1) 

x3  =  g3(c1,  c2,  c3) 

The  calibration  function  (g),  is  a  function  that  relates  the  measured 

intensity  (x)  to  the  concentration  of  analyte  (c)  via  the 

slope  of  the  calibration  curve.      The  contributions  to 

the  total  measured  intensity  at  any  one  wavelength  from  the 

various  concentrations  of  matrix  constituents,  c.  ,  c_ ,  &  c„, 

are  assumed  to  be  independent  of  one  another.   Therefore, 

x.  is  given  as  a  linear  combination  of  the  intensities  derived 

from  the  individual  concentrations  of  sample  constituents  in 

a  matrix  consisting  of  n  components.   In  general,  for  a 

multicomponent  matrix  with  n  components,  the  following  set 

of  equations  can  be  derived  expressing  the  signal  intensity, 

x.,  monitored  at  any  one  wavelength, 

xl  "  yllcl  +  ^12c2  +  : ' +  ^lncn 


X2  =  Y21C1  +  Y22C2  +  +  ^2nCn 


x   =  y  _.  c.  +  y  Oc0  +  +  Y   c 

n    ' nl  1    'n2  2  'nn  n 


(2) 
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or, 

Xj  ■   I  Y. ,C,  (3) 

3x. 
where  y   =  — -  is  the  slope  (partial  sensitivity)  of  the 

1J    3c. 
calibration  curve  when  plotting  x.  vs   c.  (x.  is  the  signal 

for  element  i  at  wavelength  A.  and  c.  is  the  concentration 

of  the  non-analyte  species  j).   It  should  be  emphasized  that 

this  equation  is  only  applicable  to  the  case  in  which  each 

component  (c^,  c2,  c3,  c  )  independently  contributes  to 

the  signal  at  the  analytical  line  (A.),  and  in  which  no 

multiplicative  correlations  appear  amongst  the  sample 

components,  i.e.,  no  intercomponent  effects  like  phosphorous 

on  calcium  (32)  in  flame  spectrometry  occur.   Actually, 

expression  (2)  can  be  applied  to  atomic  spectrometry  only 

when  the  effects  of  concomitant  species  are  limited  to 

spectral  interferences  which  are  caused  by  overlapping 

atomic     profiles  (of  emission  or  absorption)  of  the 

co-existing  species  at  the  wavelengths  of  the  analytes. 

From  equation  (2),  a  "calibration  matrix"  can  be  extracted 

and  takes  on  the  following  form, 


Yll'  Y12'  '  Yln 

Y21'  Y22'  '  T2n 


Ynl'  Yn2'  '  Ynn 


(4) 


With  the  calibration  matrix  formulated,  Kaiser  has  defined 
the  selectivity,  5,  as  follows, 


Y  •  • 


E  =  Min     n 

i=l.  ..n  S   I Y  ■•.  I  -  Iy-- 
k=l   lk     ,T±1 


-  1 


(5) 
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where  JVIin     means  that  after  evaluating  the  selectivity 

i=l .  . .  n 
with  respect  to  all  of  the  components,  1=1... n,  the  final 

value  chosen  to  be  the  indicative  selectivity  of  the 

technique  is  that  value  which  represents  the  worst,  or 

minimum,  value.   It  is  therefore  assumed  that  the  entire 

procedure's  selectivity  is  limited  by  its  worst  individual 

measure  of  selectivity.   Again,  it  should  be  stressed  that 

5  is  independant  of  the  concentration  of  the  species  'i'  in 

the  sample  solution  as  well  as  any  intercomponent  effects. 

In  analytical  spectrometry,  particularly  in  atomic 

spectrometry,  chemical  and/or  physical  interferences  (i.e., 

enhancement  or  depression  in  the  production  of  analyte  atoms 

or  ions  in  the  gaseous  environment,  such  as  flame,  plasma, 

furnace,  etc.  (.33))  are  often  more  serious  than  spectral 

interferences.   In  this  case, the  selectivity  of  an  analytical 

technique  becomes  concentration  dependent  upon  the  individual 

components  in  the  sample.   Therefore,  a  method  for  discerning 

the  selectivity  of  an  analytical  technique  when  considering 

intercomponent  effects  as  well  as  spectral  interferences  when 

either  concentration  dependent  or  independent  needs  to  be 

developed.   The  following  expression  shows  the  relationship 

between  the  analytical  signal  measured  and  any  spectral 

and  intercomponent  interferences, 

j=n         j=n 
x     =  S  y   c   +  (  S  y   c  )  (6) 

j=l  1J  J     j=l  1J  J  ci 
9x. 
where  u   =C — -).   is  the  slope  of  the  calibration  curve  of 
1J    9c.   i 
J 
x±   vs  c.  but  evaluated  at  the  concentration  of  analyte  species 
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(i)  equal  to  c.   Experimentally,  u. .  can  be  obtained  by 


observing  the  dependence  of  signal  x.  at  analyte  concentration 
c,  on  the  concentration  of  interferent  c.,  where  the 
calibration  solution  also  contains  constant  concentrations 
of  other  concomitant  species  (.1,  2,....n)  as  well.   In  this 
case,  it  is  obvious  that  all  diagonal  elements  will  by  defini- 
tion have  a  value  of  zero.   Furthermore,  equation  (6)  will 
revert  back  to  equation  (3)  if  no  inter component  effects  occur. 

In  a  similar  manner  to  equation  (2), expression  (6)  may 
be  rewritten  as  follows, 

Xl  =    Y11C1    +  (Y12+^12)C2  +  +  (Yln+Vln)cn 

X2  =  (Y21+1J21)cl  +     Y22C2     +  •'•■■  +  (Y2n+U2n)cn 

Xn  =  ^nl+ Unl)cl  +  (Yn2+V2)c2  + +    YnnV 

From  equation  (7),  just  as  we  did  for  the  spectral  interference 

case  (refer  to  expressions  (2)and(4)),  a  calibration  matrix 

of  the  analytical  procedure  can  be  given  as, 

Yn   (Yi2+lii2)  (Ym+ym) 

(Y21+li21)  Y22  (Y2n+1J2n)  (g) 


^nl^nl}    (Yn2+W    Ynn 


Equation  (8)  is  equivalent  to  the  sum  of  the  following  two 
matrices  which  represent  a  spectral  interference  calibration 
matrix  plus  an  intercomponent  interference  calibration  matrix, 

Yll   Y12  Yln  °   v12  Hn 

Y21   T22 Y2n  y21    °   y2n,n, 


Ynl   Yn2  ' Ynn  *nl   Wn2  ° 
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Having  introduced  the  intercomponent  interference  into 
our  calibration  matrix,  we  are  now  in  the  position  to  refine 
Kaiser's  definition  of  the  selectivity  of  an  analytical 
method  to  include  these  interferences  as  well  as  to  introduce 
a  concentration  dependence  upon  these  interferences.   A 
concentration  dependent  total  selectivity,  A  ,  which  includes 
both  the  spectral  and  intercomponent  interferences,  as  well 
as  the  effect  of  concentration  on  the  signal  levels,  can  be 
defined  as 

:!:y.-|c. 

Ac  =  Min    5 ^ _  i  (10) 

.  1I  'ik   Hij '  j    '  rii'  i 

Expression  (10)  accounts  for  the  effect  of  component  concen- 
trations upon  both  spectral  and  intercomponent  interferences. 
If  intercomponent  interferences  are  absent,  then  equation 
(10)  becomes  transformed  to  the  following, 

Iy11|o1 

Hc  =  Min    - 1  (H) 

i=l. . n  „   |    |      I    I 

and  if  only  the  intercomponent  interferences  are  present, 
then  equation  (10)  reverts  to 

I Y • • I c- 

?c  =  Mi?    n X  (12) 

1=1"nz   |u..|c. 

Similarly,  if  the  interferences,  both  spectral  and 
intercomponent,  are  observed  to  be  concentration  independent 
then  expressions  (10)  -  (12)  will  take  the  form  of  the 
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following  equations 


(13) 


A  ~   Min    -  1 

i=1"n  r  |y..  +  y.  ,|  -  |y.,| 
It.. I 

H  =  Min 1  (14) 

i=1--n  e   |y,,|  -  |y..| 
.  1I  'ij  '    "n1 

!y--  I 

"n1 

C  =  Min    -— 1  (15) 

i=l. .n  -   i    i 

1      'yil' 

respectively. 

In  all  cases,  £>£  ,  A>A  ,  H>H  .   For  example,  when 

Cy  C  O 

analyzing  sea  water,  the  high  concentration  of  sodium 
appreciably  interferes  with  the  atomic  absorption  signals 
of  several  elements;  in  such  a  case,  2  can  be  appreciable 
even  though  5  may  indicate  considerable  spectral  interferences, 
Finally,  it  should  be  stressed  that  for  a  completely  selective 
analytical  procedure,  H,  a  ,  A,  A  ,  %    and  £   are  all  infinite, 
whereas,  for  a  completely  non-selective  procedure,  the  same 
parameters  will  all  be  -1.   Any  finite  value  less  than 
infinity  and  greater  than  -1  will  indicate  an  analytical 
procedure  with  some  degree  of  selectivity. 


Experimental 

Materials .   All  reagents  used  were  of  reagent  grade 
quality  and  standard  stock  solutions  were  prepared  according 
to  the  procedure  of  Smith  and  Parsons  (34).   Dilutions  of  the 
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stock  solutions  were  made  covering  the  concentrations  of 
elements  as  represented  in  Table  III. 

Instrumental.   A  Perkin-Elmer  atomic  absorption  spectro- 
photometer (model  303)  was  used  for  the  flame  atomic 
absorption  spectrometry  studies.   The  air/acetylene  flame 
(C2H   flow  rate:  3-6  1/min,  air  flow  rate:  10  1/min)  was 
ignited  using  a  10  cm  slot  burner.   The  experimental  conditions, 
such  as  the  burner  position  and  fuel  flow  rate,  were  adjusted 
to  obtain  suitable  signals.   For  the  atomic  emission  plasma 
measurements,  a  SMI  Spectrospan  II  emission  spectrophotometer 
with  a  SMI  Spectrojet  III  plasma  excitation  source  was  used. 
The  detailed  conditions  of  atomic  emission  observation  were 
similar  to  references  13  and#14.   The  spectral  bandpass  for 

the  atomic  absorption  and  emission  spectrometers  were  6.8  A 

*- 
and  0.1  A  respectively. 

Procedure.   The  measured  analytical  lines  of  elements 
studied  are  listed  in  Table  III.   In  most  cases,  the  standard 
solution  was  not  contaminated  with  the  other  elements  (in 
the  concentration  range  of  interest)  except  for  the  following: 
zinc  in  potassium  chloride,  calcium  carbonate  and  magnesium 
metal  powder;  manganese  and  iron  in  magnesium  powder;  and 
copper  in  potassium  chloride  and  magnesium  powder.   The 
extents  of  these  contaminations  were  estimated  from  emission 
spectrometry,  namely,  the  estimated  emission  signals  were 
monitored  by  scanning  the  echelle  monochromator ,  and  the 
observed  emission  signals  were  measured  at  just  shorter  and 


Table  III 


Composition  of  Artificial  Sea  Water 
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Element 

K 

Na 

Mg 

Ca 

Sr 
Cu 
Zn 
Fe 

Mn 

Cd 


Concentration 
(ppm) 

Wavelength 

380 

4044 

10,500 

5890a 

1,350 

2796, 
2852° 

400 

4228 

8.0 

4607 

0.5C 

3247 

0.5C 

2139 

0.5C 

2483, 
3581° 

0.5C 

2794. 
4033° 

0.5C 

2288 

a  -  diluted  _1,000  times 

b  -  only  for  emission  spectrometry 

c  -  artificially  added 

Note:  Concentrations  express  those  of  cationic  species  only 
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longer  wavelengths  than  the  analyte  wavelength.   When  no 
emission  signals  occured  at  these  wavelengths,  the  emission 
at  the  analyte  wavelength  was  assumed  to  be  due  to  contami- 
nation, and  so  the  results  were  corrected  accordingly. 

Results  and  Discussion 

Spectral  selectivity.   Calibration  matrices  according 
to  equation  (4)  have  been  obtained  with  respect  to  atomic 
absorption  spectrometry  (AAS)  and  atomic  emission  spectro- 
metry (AES)  for  a  simulated  sea  water  analysis.   The  cali- 
bration matrices  for  the  sea  water  analysis  are  shown  in 
Tables  IV  and  V. 

For  AAS,  major  elements  such  as  K,  Na,  Mg,  and  Ca 
produce  spectral  interferences  at  analytical  lines  of  the 
other  elements  (at  trace  and  lower  concentrations).   The  main 
cause  of  these  interferences  seems  to  be  the  molecular  absorp- 
tion of  these  major  elements  such  as  the  broad  band  absorption 
of  sodium  chloride  in  the  UV  region  (35,  36,  37).   The  absorp- 
tion of  sodium  chloride  overlaps  with  the  analytical  lines 
of  several  trace  elements  (Fe,  Zn,  Cd).   Although  the  partial 
selectivity  (-y i . )  is  much  smaller  than  y±±   of  such  trace 
elements,  a  spectral  interference  due  to  10,000  ppm  sodium 
can  not  be  ignored  in  the  determination  of  sub-ppm  concentra- 
tion of  trace  elements  in  the  case  of  flame  AAS  (refer  to 
5ci  value  in  Table  IV).   Actually,  EQ±   of  AAS  for  sea  water 
analysis  is  mainly  limited  by  the  interference  of  sodium  on 
the  iron  absorption  line. 
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Surprisingly,  AES  with,  the  DC  plasma-echelle  monochromator 
system  shows  higher  selectivity  than  AAS.   For  the  manganese 
(4033  A)  and  calcium  (2288  A)  lines,  10,000  -  20,000  ppm 
sodium  gave  small  negative  signals.   This  interference  might 
be  interpreted  as  being  indicative  of  a  reduction  of  the  back- 
ground emission  of  the  plasma  due  to  the  high  concentration 
of  sodium  chloride.   Small  interferences  due  to  the  sodium 
chloride  on  Fe  and  Cd  were  observed. 

The  selectivities  calculated  from  the  calibration 
matricies  are  summarized  in  Table  VI..   The  concentration 
independent  selectivities  (S)  are  30  -  10^  times  larger  than 
the  concentration  dependent  selectivities  (E  ),  and  selec- 
tivities  for  AES  with  a  DC  plasma-echelle  monochromator 
system  are  10  -  104  times  larger  than  AAS.   This  can,  in 
part,  be  attributed  to  the  higher  resolution  of  the  spectro- 
meter used  in  AES  (approximately  10X  higher)  than  the  one 
used  in  AAS. 

Inter component  selectivity.   The  calibration  matrix  for 
intercomponent  interferences  in  the  sea  water  analysis  by 
AAS  is  shown  in  Table  VII.  When  determining  the  matrix,  the 
blank  solution,  which  includes  all  species  except  analyte, 
was  used  to  obtain  the  base  line  in  order  to  exclude  spectral 
interferences  and  contamination. 

The  major  interelement  interference  in  the  matrix  is 
due  to  sodium  chloride,  which  enhances  potassium  atomic 
absorption,  whereas  the  addition  of  sodium  chloride  suppresses 
magnesium,  calcium,  and  strontium  atomic  absorptions.   The 
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Table  yi 


Selectivity  y-alues  for  Sea  Water  Analysis 


Atomic  Absorption        Atomic  Emission 
Spectrometry  Spectrometry 

67.8  6.66  X  105 


Hc  2.94  32.3 

g  22.1 

ln  0,62 
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inter component  selectivity  factors  (E • aM£«« )  for  potassium 
(enhancement)  and  for  magnesium  C suppression)  are  both 
limited  by  sodium  in  sea  water  analysis.   The  calculated 
intercomponent  selectivity  values  are  shown  in  Table  VII. 
It  is  clear  that  the  intercomponent  selectivity  is  much 
smaller  than  the  spectral  selectivity  in  the  case  of  AAS . 

Conclusion 
More  important  than  any  specific  conclusions  that  can 
be  drawn  from  the  examination  of  the  results  of  these  experi- 
ments  is  the  fact  that  numerical  values  of  selectivities  are 
now  made  available.   This  capability  now  affords  a  very 
valuable  figure  of  merit  to  be  incorporated  into  the  pattern 
recognition  problem.   In  additon  to  all  the  other  quantifiable 
figures  or  merit  available  such  as  cost  of  analysis,  linear 
dynamic  range  (including  limits  of  detection),  sample  size, 
sample  type  (binary  classifiable  (17)),  accuracy,  precision, 
etc. ,  the  selectivity  should  prove  to  be  an  effective  de- 
scriptor for  the  pattern  recognition  problem.   What  remains 
to  be  done  in  the  future  is  an  inordinate  amount  of  data 
collection  and  reduction  before  the  marriage  of  an  analytical 
problem  and  an  appropriate  technique  can  be  realized  from 
purely  mathematical  considerations. 


CHAPTER  THREE 
THEORETICAL  LIMITS  OF  DETECTION 

Introduction 

The  groundwork  established  in  the  previous  chapter  is 
such  that  it  can  eventually  succeed  in  bringing  together  any 
analytical  problem  with  any  available  analytical  technique. 
This  capability  would  prove  to  be  an  extremely  powerful 
diagnostic  aid  for  the  analytical  chemist  when  faced  with 
the  task  of  solving  difficult  analytical  problems.   Another 
frequently  encountered  problem  faced  by  the  analytical  re- 
searcher is  one  that  concerns  itself  with  the  improvement 
of  existing  techniques  such  that  they  can  be  made  amenable 
to  the  analysis  of  smaller  and  smaller  quantities  of  analyte, 
i.e.,  improving  detection  limits.   Some  of  the  questions 
the  analyst  may  ask  while  trying  to  achieve  better  detection 
limits  are: 

i.   which  parameters  whould  I  try  to  vary? 

ii.   in  which  direction  should  the  variation  be  made 
in  order  for  me  to  realize  an  improvement  in  the  detection 
limit? 

iii.   how  much  of  an  improvement  can  I  expect  in  con- 
junction with  a  particular  finite  change  in  operating 
conditions? 
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iv.   [finally,  at  what  point  can  I  expect  not  to  gain 
any  more  enhancement  of  detection  limits,  i.e.,  what  is  the 
theoretical  limit  of  detection  that  can  be  attained  under 
the  best  of  conditions? 

There  are  many  empirical  experimental  approaches  that 
can  be  utilized  effectively  to  improve  detection  limits, 
although  they  make  no  effort  to  answer  any  of  the  afore- 
mentioned questions.   From  a  practical  viewpoint,  these 
approaches  are  based  solely  on  the  operator's  intuition 
and  past  experience,  which  may  or  may  not  be  the  most 
economical  or  efficient  approach.   If  the  analyst  "could  be 
guided  by  theoretical  considerations  that  could  quickly 
expose  the  effects  of  altering  instrumental  parameters, 
a  more  objective  approach,  to  enhancing  detection  limits 
and  understanding  the  inherent  limitations  of  the  tech- 
niques can  be  realized.   Although  this  type  of  information 
can  be  extracted  from  the  following  theoretical  calcula- 
tions of  limits  of  detection,  this  was  not  the  intended 
purpose  of  this  work.   A  comprehensive  review  of  the 
application  of  signal-to-noise  theory  for  these  purposes 
can  be  found  in  Alkemade  et  al.  (.38)  and  Boutillier  et  al. 
(391.   Primarily,  the  objective  of  this  work  was  to  make 
a  comparison  between  the  theoretical  limits  of  detection 
of  a  variety  of  atomic  plasma  spectrometric  techniques, 
and  to  assess  the  current  ''state  of  the  art"  with  respect 
to  realizing  these  limits.   This  will  be  accomplished  by 
considering  the  major  contributions  and  factors  affecting 
both  the  signal  production  from  the  analyte,  and  limiting 
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noise  sources  associated  with  the  measurement  process. 
Furthermore,  a  theoretically  calculated  limit  of  detection 
may  prove  to  be  an  important  figure  of  merit  to  be  used 
in  the  pattern  recognition  problem  of  Chapter  Two.   In 
effect,  if  no  state  of  the  art  technique  is  available  with 
sufficiently  low  detection  limits,  the  analyst  may  choose 
to  spend  time  optimizing  if  it  is  theoretically  possible 
to  achieve  the  appropriate  sensitivity. 

To  be  more  specific,  by  means  of  signal-to-noise 
expressions  for  atomic  absorption  flame/plasma  spectrometry 
(with  narrow  line  and  continuum  sources),  for  atomic 
emission  flame  and  ICP  spectrometry,  for  atomic  fluorescence 
flame  spectrometry  with  "narrow"  line  conventional  and  high 
intensity  (pulsed  dye  laser)  sources,  for  atomic  fluo- 
rescence flame  spectrometry  with  "continuous"  conventional 
and  high  intensity  sources  (pulsed  and  cw  dye  lasers) ,  for 
photoacoustic  spectroscopy  with  a  cw  high  intensity 

source,  and  for  optogalvanic  spectroscopy  with  a  pulsed 

3 

high  intensity  source,  detection  limits  (in  atoms/cm  of 

gases)  are  estimated.   In  all  cases,  it  is  assumed  that 
procedural  errors  are  minimized,  and  that  only  noise 
sources  in  the  instrumental  methods  are  dominant .   In 
particular,  the  limiting  noises  discussed  are  source 
transmitted  shot  noise,  emission  background  shot  noise 
and  detector  dark  current  shot  noise;  and  emission  back- 
ground shot  noise;  for  atomic  absorption,  atomic  fluorescence 
and  atomic  emission  spectrometry,  respectively.   For 
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the  analysis  of  photoacoustic  spectroscopy,  Brownian  motion 
of  ambient  air  molecules  is  the  limiting  nois.e  source, 
whereas  shot  noise  due  to  the  finite  electron  population  in 
the  flame   is  taken  to  be  the  limiting  noise  source  for 
optogalvanic  spectroscopy. 

Conversion  of  the  detection  limits  in  terms  of  analyte 
atomic  (.ionic)  species/cm^  of  hot  gases  to  analyte  solution 
concentration  aspirated  into  the  flame  or  plasma,  will  be 
performed.   The  expressions,  and  particularly  the  gaseous 
detection  limits,  are  discussed  in  terms  of  the  present 
"state  of  the  art'1'  of  trace  analysis,  and  with  respect  to 
expected  improvements  in  detection  limits  as  well.   Also, 
the  conversion  to  detection  limits  in  liquid  solution  is 
discussed  particularly  with  respect  to  the  agreement  between 
the  predicted  detection  limits  via  the  signal-to-noise 
ratio  expressions,  and  the  experimentally  measured  detection 
limits. 

Explicit  derivations  of  the  signal  and  noise  expres- 
sions will  not  be  presented  here,  although,  for  in  depth 
treatments, the  reader  is  referred  to  the  appropriate  ref- 
erences as  they  are  cited  in  the  text.   Furthermore,  due  to 
the  natural  division  amongst  the  techniques  that  will  be 
analyzed,  i.e.,  those  that  use  photomult iplier  tubes  as  a 
detection  device  and  those  that  do  not  (microphone  and 
electrodes),  the  former  techniques  will  be  treated  sepa- 
rately, holding  off  the  analysis  of  photoacoustic  spectros- 
copy and  optogalvanic  spectroscopy  until  the  end. 
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The  detection  limit  will  be  defined  according  to  the 
IUPAC  recommendation  (29),  i.e.,  the  concentration  resulting 
in  a  S/N  =  3.   For  all  atomic  spectroscopic  techniques,  S/N 
expressions  will  be  given  in  terms  of  digital  integrated 
values  (.counts);  all  "intensities"  will  be  given  in  terms 
of  counts  per  second  per  x  per  y  per  z  (x,  y,  z  depend 
upon  the  type  of  intensity:  an  irradiance  is  in  counts/s  cm2, 
a  spectral  irradiance  is  in  counts/s  cm  nm,  a  radiance  is 
in  counts/s  cm2  sr,  a  spectral  radiance  is  in  counts/s  cm2 
sr  nm,  etc.).   The  S/N  expressions  will  be  given  so  that 
the  numerator  is  the  complete  signal  expression  and  the 
denominator  is  the  complete  limiting  noise  expression. 

In  the  following  analysis,:  only  atomic  spectral 
techniques  involving  plasmas  (flame,  inductively  coupled 
plasmas,  etc.)  will  be  considered.   No  attempt  will  be  made 
to  consider  detection  limits  in  furnaces,  arcs  and  other 
devices  where  solid  samples  are  measured.   The  reasons  for 
this  differentiation  are.  i)  the  atomic  plasma  spectrometric 
methods  all  involve  nebulizers  and  production  of  gaseous 
species  in  the  plasma,  ii)  for  solid  sample  devices,  esti- 
mation of  the _ concentration  of  analyte/cm3  is  of  little 
interest,  but  rather  the  estimation  of  the  absolute 
amount  (in  g)  of  the  analyte  to  produce  a  certain  species 
concentration  is  of  greater  interest,  iii)  solid  sample 
devices  produce  a  transient  change  in  analyte  species,  and 
so  both  the  integrated  signal  (.over  time)  and  the  peak 
signal  (at  the  time  of  maximum  response)  should  be 
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calculated,  iv)  the  noise  sources  in  furnaces  are  not  as  well 
documented  as  flames  and  plasmas.   In  conclusion,  due  to  the 
aforementioned  points,  there  is  no  fair  means  of  comparing 
furnace-arc  atomic  spectrometric  systems  with  plasma  atomic 
spectrometric  systems,  and  so  only  the  latter  atomic  systems 
will  be  considered  here.   We  are  now  prepared  to  discuss  the 
S/N  expressions  in  more  detail. 

Analysis  of  Photomultiplier  Techniques 
Atomic  Fluorescence  Flame  Spectrometry  (AFFS) 

It  will  be  assumed  here  that  the  major  (limiting) 
noise  in  AFFS  is  either  detector  shot  noise  (an  ideal 
situation)  or  flame  background  emission  shot  noise,  i.e., 
such  noises  as  detector  flicker  noise,  background  emission 
flicker  noise,  source  scatter  shot  and  flicker  noises, 
band  (or  line)  concomitant  fluorescence  shot  and  flicker 
noises,  band  (or  line)  concomitant  emission  shot  and 
flicker,  and  electronics  noise  are  minimized  either  by 
modulation  techniques,  by  physical  separation  of  the 
concomitants,  by  more  efficient  atomization  procedures,  or 
by  specialized  methods  such  as  ratioing  techniques.   The 
detector  shot  noise  or  the  background  emission  shot  noise 
are  the  remaining  irreducible  and  limiting  noises  under 
these  conditions. 
Continuum  source  excitation  (AFFSC) 

The  limiting  detectable  S/N  for  AFFSC  assuming  a  low 
intensity  continuum  source,  CL,  is  given  by  (40): 
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A„        i?         fy      ^     xr     a  U        L, 


AFESCL     10'  D   1  %      E   (X  )  Y  A   -£  n^UwhT  (Le-i,.):. 

S  4.7T  S^  PA   °      Ul  gl   T      s  D  D  E   i 


S/N) 

lim,  DS 


C2  DD  x.  rd  nE)* 


assuming  detector  siiot  noise,  DS ,  is  limiting  and  by 

4  2" 

AEESCL     107  D  I      Xo      Em  (A  )  Y  A  ,  -H  nL(ftwhT  Qn£nnw)x. 
S/N)         =       s  4u  87c"  PX   °      ul  gl   T      s  D  D  E   x 
lim,  BS __ 

(2  DB  BPBX^o^RwhTsVDT1BKAA8)Ti>* 
assuming  background  shot  noise,  BS ,  is  limiting.   All  terms 
are  defined  below  : 

Aul  =  transition  probability  (Einstein  coefficient  of 

spontaneous  emission)  for  emission  process,   s 
I        =  emission  (fluorescence)  path  length,   cm 
XQ      =   peak  absorption  wavelength,   cm 
Y   =  quantum  efficiency  of  transition,  dimensionless 
g's  =  statistical  weight  of  upper,  u,  and  lower,  1,  levels, 
dimensionless 

n   =  limiting  detectable  concentration  of  analyte  in 

T  -3 

flame,  cm 

Q   =  spectrometer  or  entrance  optics  solid  angle 

(whichever  is  limiting) ,   sr 
w,h  =  width,  height  of  spectrometer  slit  (assumed  to  be 

limiting  area),  cm,  cm 
Ts  =  spectrometer/entrance  optics  overall  transmission, 

dimensionless 

s   =  duty  factor  for  excitation/measurement  of 
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fluorescence,  dimensionless 
DD  =  duty  factor  for  detector,  dimensionless 

Ep^CAQ)  =  spectral  photon  energy  at  central  absorption 

—2   -1    —1 
wavelength,  cm   s   nm 

DB  =  duty  factor  for  measurement  of  Background  emission 

in  flame,  dimensionless 
QD  =  photocathode  quantum  efficiency  of  photomultiplier 

detector,  dimensionless 
£D  =  first  dynode  efficiency  of  photomultiplier  detector, 

dimensionless 
n.E  =  efficiency  of  amplifier  -  discriminator  -  readout 

system,  dimensionless 
Rq  =  average  dark  current  rate  of  detector,  counts  s 

BpB^(.Ao)=  average  flame  emission  background  photon  radiance, 

-1    -2    -1   -1 
counts  s   cm   sr   nm 

AA  =  spectrometer  bandpass,  nm. 
The  limiting  detectable  S/N  for  AFFSC  assuming  a  high 
intensity  (CH)  continuum  source  (I.e. ,  one  producing  satura- 
tion) is  given  by, 

AFFSCH     D  I      A    (fTF")  nT  (»T  Q-S-nJi, 
S/N)         =    s   4ir      ul      gu  gl    T       s  D  D  e   i 


lim,  DS 


<2  DD  Ti  *D  V 
assuming  detector  shot  noise,  DS ,  is  limiting  and  by, 


AFFSCH 
S/N) 

lim,    BS 


gu    N      L 


.    Ds  I  \i  CiTi^  %  CP-WA^n,)^ 


(2   DB   AAs   BPBX<VCnwhTs(SDeDVTl>2 
assuming  background  emission   shot   noise,    BS ,    is   limiting, 


58 


Line  source  excitation 

The  limiting  detectable  S/N  for  AFFSL  assuming  a  low 

intensity  narrow  line  source  (_LL)  ,  i.e.,  one  with,  a  line 

wdith,  much,  narrower  than  the  absorption  line  width  is-  given 

by, 

4  g 

AFFSLL        107D      £     %__  BpCA>T  A        JI  nj  £,    CSwM  %^^)i. 
S/N)  =  s    4tt    8.tt.c      P      °  ui    gl      T    6XD  s   D  D   E      1 


lim,    DS 


2 


(2  DD  Ti  *D  V 
assuming  detector  shot  noise  is  limiting  and  by, 


,7,  4 


AFFSLL        10    D       £     Xo      ILU    )Y   A   -    -£  n£  u      (fiwhT   Q-e-.n-K. 
S/N)                  =  s   4tt  Bire     P      °  "        Ul    gl      T    6*'  s   D  D   E      x 

lim,    BS      D 


<2DB  ^PB^o^^sVdV**   V 


on 


assuming  background  emission  shot  noise  is  limiting.   All 
terms  above  have  been  previously  defined  except  for, 

EpCVQ)  =  line  source  photon  irradiance  at  the  absorpti 

-1    -2 

line  center,  X  ,  photons  s   cm 

u  ■  profile  factor  (.41)  accounting  for  relative 
width  of  source  line  and  absorption  line  and 
for  their  profiles,  dimensionless 
When  using  a  high  intensity  (.laser)  narrow  line  source,  two 
general  cases  can  be  considered,  namely  interaction  of  the 
laser  line  with  a  homogeneously  broadened  absorption  line 
(collisional,  Lorentzian  broadening)  and  interaction  of  the 
laser  line  with  a  small  group  of  absorbers  having  the 
proper  translational  (Doppler  interval)  energy  which  is  the 
case  of  inhomogeneous  broadening.   The  S/N  expression  for  the 
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latter   case   is  unique   and   is   given  below.      The   limiting 

detectable   S/N   for   AFFSL   assuming   a  high,  intensity   line 

source,    LH?    is   given  by, 

T.  XQ7E^{X    )X   4   , 

AFFSLH        D      I      A        n^C P      °      °    )3(flwiiT   Q   e    n.J.T, 

S/N)                  =      s   4tt      u1       r   16ttc    SX     ..  s   D  D   E      i 

lim,    DS     ett 


C2DD   Ti   RD   V 

assuming  detector  shot  noise  is  limiting  and  by, 

T  10  V, .(X  )X  4  . 

AFFSLH   D   1  A   n£( 2-2-i.)t(fllrtlT  o  e  n  )t 

S/N)        =   s  4tt   Ul   i  i67TC  5X             s  D  D  E   i 
lim,  BS  ™ 

(-2DB  BPBX<V^whTs%£DVAAs  Ti)4 
assuming  background  emission  shot  noise  is  limiting.   All 
terms  have  been  defined  previously  except  for: 

^eff  =  ^^ective  absorption  linewidth  =  tt/2  5Xt  ,  nm 
<5XL   =  Lorentzian  absorption  line  width,  nm 
The  preceeding  expression  for  atomic  fluorescence  flame 
spectrometry  are  specifically  meant  to  apply  to  the  case 
where  the  typical  atom  considered  is  a  two-level  one  under- 
going a  single  photon  excitation.   However,  two-photon 
processes  are  showing  significant  analytical  utility  (42,43) 
and  hence  will  be  considered  here  as  well.   The  two  processes 
of  major  importance  are  the  two -photon  absorption  via  a 
virtual  intermediate  level  and  that  via  a  real  intermediate 
level.   For  the  case  of  the  virtual  level  the  S/N  expression 
is  given  by  the  following: 
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2   photon,    virtual      Ax      I      A        n^  E2    ( A    )  S   YCQwhT  Q^twJt  . 
S/N)  L  Ttt  a  D  D  E     i 


11m,    DS 


assuming  the  use  of  a  low  intensity  continuum  source  of 
excitation  and  detector  shot  noise  limitation.   All  terms 

have  been  previously  defined  except  for   5  ■  two  photon 

4 
absorptivity  (cms)  and  ArT '  =■ pulse  width  of  the  laser  (s). 

For  the  case  in  which  a  real  intermediate  level  is  used,  the 

S/N  expression  is  given  by  the  following  (assuming  g  »  g„ ) , 

S/N)'  Ph0t°n'  "^^   D*  ^  ^  ^   (^WhTsQD^E)Ti 


lim,  DS 


C2DD  Ti  RD  V 


where  all  terms  have  been  previously  defined  and  assuming 

a  high,  intensity  continuum  source  of  excitation  that  produces 

saturation  of  both  of  the  two  photon  transitions. 

Atomic  Emission  Plasma  Spectrometry  (AEPS) 

In  this  case,  the  plasma  can  be  a  dc  plasma,  an 
inductively  coupled  plasma,  a  capacitively  coupled  microwave 
plasma,  a  flame,  etc.    The  major  noise  source  will  be  assumed 
to  be  background  emission  shot  noise,  although  in  reality 
background  flicker  noise  or  emission  shot  noise  or  flicker 
noises  due  to  concomitants  may  be  the  major  noise  source. 

Assuming  the  background  shot  noise  limitation,  the 
limiting  S/N  expression  is  given  by  the  following, 
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I    4  g 
S/N)_AEPS   „  10?  l   8^  g  V  »T  gX^WB^ 

lim,  BS  —       ~  ~ 

l2BPBA(Xo)('QwhTsQD^E)AXs  Ti 

where  all  terms  Lave  Been  defined  previously  except  for 
the  following, 

BPBA^o^  ~  averaSe  photon  spectral  radiance  of  plasma 

Background  evaluated  at  the  analyte  line 

-1    -2    -1    -1 
center,  photons  s   cm   nm   sr 

Bp^(.^0)   =  average  photon  spectral  radiance  of  a 

blackbody  radiator.,  evaluated  at  the  analyte 

line  center  and  at  the  plasma  temperature, 

photons  s~   cm   nm   sr 

I  =   emission  pathlength,  cm. 


Atomic  Absorption  Flame  Spectrometry  (AAFS) 

In  atomic  absorption  flame  spectrometry,  the  major  noise 
sources  are  associated  with  source  transmitted  shot  and 
flicker  noise  and  flame  background  emission  shot  noise. 
Other  noise  sources  can  generally  be  minimzed  by  various  types 
of  modulation  techniques  or  are  negligible  compared  to  the 
above  major  noise  sources.   In  this  analysis,  we  will  assume 
the  major  noise  source  is  source  transmitted  shot  noise,  SS . 

Continuum  source  excitation  (AFFSC) 

In  this  case,  the  limiting  S/N  expression  is  given 
by  the  following  equation, 
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AAFSC        _   107Ds    1  ^  !a  Aul   n£  ^HMI,^,),, 


11m,      SS 


C2Ds  WV'^sVdV^b  Ti^ 


Line  source  excitation  (AAFS'L) 

For  the  case  of  line  source  excitation  the  S/N 

expression  is  given  by, 

,  4 

S/N)AAFSL  .  10\   '4|\!  °T^K,C0^r.^«Bn,)T1 

lim,  SS  r 

(2Ds  BpLCOwLTg.QDeDnE)ti).^ 

where  all  terms  except  those  below  have  been  previously 
defined, 

BPL  =  Photon  radiance  of  line  source  evaluated  at 

1     —2     1 
analyte  line  center,  photons  s   cm   sr 

Evaluation  of  Other  Noise  Sources 

In  the  above  cases  only  specific  noise  sources  for 

each  method  were  considered  and  the  specific  noises  (,in 

counts)  were  given  in  the  denominator  of  each  expression. 

In  this  section,  several  other  noise  expressions  Cin  counts) 

are  given  and  can  be  substituted  in  the  denominator  for  the 

specific  case  if  one  wishes  to  obtain  S/N), .   for  other 

'       lim 

limiting   cases   than   those   given  previously. 

Plasma    (frame)   background   flicker  noise 

NFBF   =    <-2^    ^B  5PBAao)(^whTsQD£DnE)AAs    Ti 
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where  §_  *   flicker  factor  for  background  emission 

as  explained  in  Appendix-  B?  dimensionless 


Continuum  source  flicker  noise  in  AA 

Ncp  -  C21f  Cc  BpCAa0KnwnTS:QDsDnEl(A'Asl  x. 
where  ?c  =  flicker  factor  for  continuum  source, 
dimensionless 

Line  source  flicker  noise  in  AA 

NLF  =  C2i*  SL  BLCXo)(awhTgQDeDnE)Ti 

where  £_  =  line  source  flicker  factor,  dimensionless 

Continuum  source  scatter  flicker  noise  in  AF 

NSCF  =    C2>*    5SC    (^S}    nS   *   EPSCCVWwhTs'VDVTi 

where  =  differential  cross-section  for  scatter, 

d^       2    -1 
cm   sr 

?sc  ■  continuum  source  flicker  factor, 

dimensionless 

ns   =  concentration  of  scattering  species, 

-3 

cm 

EpSCCA0)   =  spectral  source  (.continuum)  irradiance 

-1   -2   -1 
at  X    ,  photons  s    cm   nm 

I        =  pathlength  of  scatter  in  direction  of 

detector,  cm 
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Line  source  scatter  flicker  noise  in  AF 

NSLF  *  ^    %L   EPLao)-(^s1  nS  ^:WTsQDsDr1.E).Ti 

dft 
where  ^SL  =  line  source  scatter  flicker  factor, 
dimensionless 

EpLCAQ).  -  line  source  photon  irradiance  at  \    , 

-1    -2 
photons-  s   cm 


Continuum  source   scatter   shot   noise   in   AFS 


scs  =  ^psc^s^s  K^hTsQD£DnE)AAs  x.) 


Line  source  scatter  shot  noise  in  AFS 

NSLS  "  t2EpL/.d^ns  £CfiwhTsQD£DnE)Ti)* 
dft 
No  expressions  will  be  given  for  concomitant  band  or 
line  fluorescence,  emission  flicker  or  shot  noise,  because 
these  noise  sources  are  very  specific  (to  certain  wavelength 
regions)  and  so  can  be  avoided  either  by  separation  in  a 
spectral  region  free  of  concomitant  fluorescence  or  the 
concomitant  can  be  physically  separated  and  thus  be  made 
negligible.   Finally,  no  noise  expressions  will  be  given 
for  multiplicative  noises,  such  as  variation  in  aspiration 
flow  rate,  droplet  dispersion  rate,  desolvation  variations, 
and  vaporization  changes  because  they  can  be  minimized  by 
either  instrumental  improvements  and/or  by  the  use  of 
internal  standards  or  ratioing  methods. 
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Results 

The  results  of  th.e  calculations,  of  theoretical  detection 
limits  for  various  plasma  spectrometric  techniques-,  employing 
photomultiplier  tubes-  as  the.  detection  device,  can  he  found 
in  Tables  VI IT  and  IX.   The  parameters  used  to  estimate  the 
detection  limits  are  listed  in  Table  X.   The  values  of  the 
parameters  were  chosen  in  the  following  manner: 

i)  for  all  cases,  a  highly  probable  transition  (high 
A^^r'was  chosen; 

ii)  for  all  cases,  a  short  (.300  nm)  and  a  long  (.600  nm) 
line  was  chosen  to  indicate  the  effect  of  wavelength; 

iii)  for  all  atomic  fluorescence  calculations,  realistic 
estimates  of  quantum  efficiencies  were  made  for  the 
flames/plasmas  of  interest ; 

iv>  flame  background  spectral  radiance  estimates  were 
based  upon  the  results  of  Gilbert  (.44)  and  Fu j  iwara  et  al . 
(45).   ICP  plasma  background  spectral  radiance  estimates 
were  based  upon  emissivity  estimates  from  Dresvin  (.46)  and 
temperature  measurements  by  Eckert  (47); 

v)  Doppler  half  widths  were  calculated  (48)  assuming 
an  atom  of  50  amu  and  the  temperature  of  the  flame  or 
plasma;  the  Lorentzian  half  width  was  estimated  from  the 
work  of  Wagenaar  (49);  the  effective  lifetime,  x,  was 
estimated  from  Lijnse  (.50); 

vi)  blackbody  spectral  radiances  came  from  Planck's 
radiation  distribution  law  as-  listed  by  Pivansky  and  Nagel 
(51); 
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Table  IS 


Calculated  Detection  Limits  for 
Two  Pilot  on  Excitation  of 
Atomic  Fluorescence  Flame  Spectrometry 


Virtual  Level 


Concentration 
(atoms/cm  ) 


'"'  s"1)         2.0  X  106 
Detector  Noise  Limit  (R_  =  104  s-1)         2.0  X  107 


Detector  Noise  Limit  (E     =  10   s   )         2.0  X  10 


Real  Level 

Detector  Noise  Limit  (R  =  102  s"1)         8.5 


D 
*D 


Detector  Noise  Limit  (_R  =  104  s"*1)         8.5  X  101 
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Table  X 


Evaluation  of  Parameters  used  to  Calculate  Detection  Limits 

(Species/cm3) 


Photon/Atom-Ion  Characteristics 

Aq    =   300  nm;  hvQ   =   6.6  x  10"19J 
Aq    =   600  nm;  hvQ   =   3.3  x  10_19J 

Y  =  0.01  (atoms  in  C2H2  -  based  flames) 

Y  =0.5   (atoms  in  H2/0  /Ar  flame) 

Y  =1    (atoms /ions  in  ICP) . 

8   —  ~\ 
Su        =  Sj  '<    A   =10   s    (for  all  atoms/ions  assuming 

highly  probable  transitions) 

5A  ._  =  0.001  nm 
ef  f 

5AD(Ao  =  300  nm,  M  =  50  amu)  =  0.0015  nm  (T=2500  K) 
<5AD(Ao  =  600  nm,  M  =  50  amu)  =  0.003  nm  (T-2500  K) 
5Ad(Aq  =  300  nm,  M  =  50  amu)  =  0.0016  nm  (T=3000  K) 
5Ad(Aq  =  600  nm,  M  =  50  amu)  =  0.003  nm  (T=3000  K) 
5AD(Ao  =  300  nm,  M  =  50  amu)  =  0.0024  nm  (T=6000  K) 
5Ad(Aq  =  600  nm,  M  =  50  amu)  =  0.005  nm  (T=6000  K) 
u   =  0.5 
Plasma/Flame  Characteristics 
C2H2/Air  Flame:  T=2500  K 

BBX  (300  nm)  =  10-6  J  s"1cm"2sr"1nm"1( 1 . 5  x  1012 

-1   -2   -1   -1 
photons  s   cm   sr   nm   ) 


B 


-1-2   -1   -1,«    h„10 


BX  (600  nm)  =  10  '  J  s   cm"  sr   nm   (3  x  10 

-1   -2  -1       -1 
photons  s   cm   sr   nm   ) 


continued 
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Table  X-continued 

C2H2/N20  Flame:  T=3000  K 

§BX  (300  nm)  -  10~5  J  s~1cm~2sr""1nm~1(  1 . 5  x  1013 

,  .      -1-2   -1   -1N 
photons  s   cm   sr   nm   ) 

BBA  (600  nm)  =  10~7  J  s~1cm~2sr"1nm~1( 3  x  1011 

.   .      -1   -2   -1   -1, 
photons  s   cm   sr   nm   ) 


H2/02/Ar  Flame:  T=2500  K 


BBX  (300  nm)  -  10  6  J  s  1cm  2sr  Xnm  1(1.5  x  1012 

,  .      -1-2   -1   -1. 
photons  s   cm   sr   nm   ) 

§BA  (600  nm)  -  10~10  J  s~1cm"2sr~1nm~1( 3 .  x  108 

u  4-     -1-2   -1   -1. 
photons  s   cm   sr   nm   ) 

ICP:  T=6000  K,  e.  =  10_5(emissivity  at  A   =  300  nm 

A  O 

and  at  600  nm) 

B  ,(300  nm)   =  10~5  J  s~1cm~2sr"1nm~1 ( 1 . 5  x  1013 

-1   "2„  -1   -Is 
photons  s   cm   sr   nm   ) 

1L,(600  nm)   =  10~5  J  s~1cm~2fer~1nnf1(  3.  x  1013 

DA 

.  ,      -1-2   -1   -1, 
photons  s   cm  sr   nm   ) 

Blackbody : 

T=2500  K 

Bbb(300  nm)  =  2  x  10~5J  s~1cm~2sr~1nm~1 

A   , 

13  -1   -2   -1   -1 

(3.  x  10   photons  s   cm   sr   nm   ) 

Bbb(600  nm)  =  1  x  10~2  J  s~1cm~2sr~1nm""1 

,„     1r.16  u    ,  -1  -2      -1   -1. 

(3.  x  10    photons  s   cm   sr   nm   ) 


continued 
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Table  X-continued 

T=3000  K 

Bbb(300  nm)  =  6  x  10"4  J  s"1cm"2sr"1nm~1 

(9.  x  10   photons  s~  cm"  sr~1nra""1^ 
Bbb(600  nm)  =  5  x  10~2  J  s~1cm"2sr~1nm"1 

17  —1   —2   — 1    1  "■> 

(1.5  x  10   photons  s   cm   sr   nm   ; 

T=6000  K 

Bbb(300  nm)  =  1.5  x  10°  J  s"1cm"2sr"1nm"1 

,18  ^  ,_   -1   -2,.  -1   -1 


(2.  x  10   photons  s   cm"  sr"  nm"  ) 

0  _   -1   -2   -1   -1 
J  s   cm   sr   nm 

~]  Q  — .  "i      _  O     _  1        "1 

(9  x  10    photons  s   cm   sr   nm   ) 


Bbb(600  nm)  =  3.  x  10°  J  s~1cm"2sr~1nm"1 


Cell  Factors 


Fluorescence  (all  flames  and  TCP): 

1  =  l^cm  (EDL;  xenon  arc) 

1  =  0.5  cm  (Dye  lasers  beam  expanded) 
Emission : 

1  =  2  cm  (all  flames) 

1  =  1  cm  (ICP) 
Absorption : 

L  =  10  cm  (C2H2-based  flames) 

L  -  5  cm  (H2/02/Ar  flame) 

L  =  1  cm  (ICP) 


continued 
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Table  X-continued 

Instrumental  Measurement  Factors 

Fluorescence:  OTh  T  Qn  ennT  =  5  x  10~4cm2sr: 
AA   =  1  nm  (all  sources  except  cw 


dye  laser) 


-5 


Emission : 


Absorption 


PWh  Tg  QD  £DnE  =  5  x  10   cni  sr; 
AA   =  0.1  nm  (cw  dye  laser  only) 
nE  =  1  (in  all  cases) 
P.Wh  Tg  QD  eDnE  =  1  x  10"6cm2sr; 
AA   =  0.01  nm  (flame  or  ICP) 
S7Wh  Tg  QD  eDnE  =  5  x  10~5cm2sr; 


AA   ■  0.1  nm  (HCL) 


-7 


OTh  Tg  QD  sDnE  =  1  x  10  'an  Sr; 

AA   =  0.01  nm  (xenon  arc  source) 
Integration  Time,  T.  =  10  s  (all  cases) 
Modulation-Pulsing  Conditions 

Fluorescence:  EDL  or  Xenon  Arc  Source: 


Ds  =  0.5,  DB  =  DD  =  1 
Modulated  CW  Ar  -Dye  Laser: 


Emission : 


Absorption 


Ds  "  °"5>  DB  =  DD  "  X 
Pulsed  N„-Dye  Laser: 

Ds  "  DB  =  DD  =  10"7 

Pulsed  flashlamp  Pumped  Dye  Laser: 

Ds  =  DB  =  DD  =  10"5 

Ds  =  DB  =  DD  =  1  (flames  and  ICP) 

D   =  0.5  (all  sources  and  all  cell  types) 


continued 
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Table  X-continued 
External  Sources  for  AA/AE 
HCL  (u  =  0.5) : 

BT  (300  nm)  =  10~4  J  s~1cm~2sr~1 

Li 


,-     -         iri14      .      ,  -1      -2  .    -1. 

(1.5   x   10        photons   s      cm      sr      ) 

§L   (600   nm)    =    10~4   J   s^cnf  2sr-1 

,„.        ,-14      .  -1      -2,     -1. 

(3   x   10        photons   s      cm      sr      ) 

Xenon   Arc    (300  W  EIMAC  Arc;    Q        =    3tt   sr; 

v  '   ex        ' 

E  ,  -  Q      B  .  ) 
ca     ex  cX 

B  v(300  nm)  =  2  x  10~2  J  s~1cm~2sr~1nm~ 

i  ^  TOT 

(3  x  10    photons  s   cm   sr   ) 

E  ,(300  nm)  =  2  x  10_1  J  s_1cm~2nm~ 
cA 

17  —1   -2 

(3  x  10   photons  s   cm   ) 

B  —1     -1   -2   -1   —1 

cA(600  nm)  =  1  x  10    J  s   cm  sr   nm 

17  -1   -2   -L 
(3  x  10    photons  s   cm   sr   ) 

Eca(600  nm)  =  1  x  10°  J  s_1cm"2nm_1 

18  —1   —2 
(3  x  10   photons  s   cm   ) 

EDL  (150  W;  n    -  0.1  Sr;  ET  =  Q      BT ,  u  =  0.5) 

ex  Li    ex  Li 

BL  (300  nm)  =  10~3  J  s"1cm"2Sr"1 

(1.5  x  10   photons  s~  cm"  sr   ) 

EL  (300  nm)  =  10~4  J  s_1cm~2 

,.    -    ,A14   .  .      -1   -2, 
(1.5  x  10   photons  s   cm   ) 

BT  (600  nm)  =  10"4  J  s"1cm~2sr~1 

Li 

(3  x  1014  photons  s_1cm~  sr~  ) 

EL  (600  nm)  =  10~5  J  s-1cm~2sr_1 

13  —1   —2 

(3  x  10    photons  s   cm   ) 


continued 
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Table  X-continued 

Dye  Laser,  CW  Ar  -Dye  Laser,  1  W  output  at  600  nm, 
0.05  W  output  at  300  nm,  beam  diameter  before 
expansion  »  0.2  cm,  beam  diameter  after  expansion  = 
0.5  cm,  band  pass  ^0.01  nm,  continuum  source  case). 

E  ,  (300  nm)  =  2.5  x  101  J  s~1cm~2nm~1 
cX 

19  -1   -2   -1 

(8  x  10    photons  s   cm   nm   ) 

-  2     -1   —2   —1 

E  .  (600  nm)  =  5  x  10   J  s   cm   nm 

C  A 

ro  -m20   .   .       -1-2   -1N 

(8  x  10   photons  s   cm  nm   ) 

Dye  Laser,  Pulsed  Np-Dye,  50  kW  peak  power,  5  ns  pulses, 

20  pps,  beam  diameter  before  expansion  =  0.15  mm, 
beam  diameter  after  beam  expansion  =  0.5  cm,  band  pass= 
0.03  nm,  doubling  efficiency  =  10%  continuum  source 
case) 

E  .  (300  nm)  =  1  x  106  J  s"1cm"2nm"1 
cX 

24  -1   -2   -1 

(1.5  x  10  ~  photons  s   cm   nm   ) 

E  .  (600  nm)  =  1  x  10   J  s-1cm~  nm" 

(3  x  10    photons  s   cm  nm   ) 

Dye  Laser,  Pulsed  N„-Dye,  same  as  above  but  with  etalon 

narrowing,  band  pass  is  now  0.001  nm,  narrow  line 
source  case) 

EL  (300  nm)  =  7  x  103  J  s-1cm~2 

22  —1   -2 

(1  x  10   photons  s   cm   ) 

EL  (600  nm)  =  7  x  104  J  s_1cm-2 

(1  x  10""   photons  s   cm   ) 

Dye  Laser,  Flashlamp  Pumped,  5  kW  fundamental,  1  ys 
pulses,  10  pps,  beam  diameter  before  expansion  -0.2cm, 
beam  diameter  after  expansion  -0.5  cm,  band  pass  « 
0.1  nm  doubling  efficiency  =  10%  continuum  source  case) 

E  ,  (300  nm)  =  1  x  104  J  s"1cm'"2nm~1 

cA  v       ' 

22  -1   -2   -1 

(1.5  x  10    photons  s   cm   nm   ) 


79 


Table  X-continued 

-  5-1-2   -1 

E  ,  (600  nm)  ■  1  x  10   J  s   cm   nm 


'X 


23  -1   -2   -1 

( 3  x  10    photons  s   cm  nm   ) 


Saturation  Conditions 

Continuum  Source 

V  *  /Qnn    %  s.  6  x  104  ,   -1   -2   -1 
E  ,  (300  nm)  >  = J  s   cm  nm 

ca        *  Y 

99 

-      ,    9  x  10zz   .  .      -1-2   -1, 
(  = —   photons  s   cm  nm   ) 

5  *  ^rtn    v  ^  2  x  103  T   -1   -2   -1 
E  .  (600  nm)  >  = J  s   cm  nm 

CA  I 

,  6  x  1021  ,  ,      -1-2   -1. 
(  y photons  s   cm  nm   ) 

Line  Source  (assuming  AXT  ^0.003  nm) 

EL*  (300  nm)  >  ~   J  s~1cm"2 

19 
.    4   x   10  .     .  -1      -2, 

(   y photons   s      cm      ) 

EL*    (600   nm)    >   i~  J   s_1cnf2 

,  2  x  1018  .  ,      -1   -2, 
(  ^ photons  s   cm   ) 
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vii>  spectral  radiances  of  the  xenon  arc  came  from 
Gilbert  (52);  line  radiances  of  electrodeless  discharge 
lamps  came  from  Mansfield   et  al.  (53);  line  radiances  of 
hollow  cathodes  came  from  Prugger  et  a'l .  (54); 

villi  all  spectral  radiances  and  line  radiances  of 
the  three  lasers  (N2-  pumped  dye  laser,  flashlamp  pumped 
dye  laser,  and  cw  Ar+  pumped  dye  laser)  came  from  manu- 
facturer's estimates  C55,  56,  57)  of  (peak)  power  and 
spectral  bandwidth;  all  laser  beams  are  assumed  to  be  ex- 
panded to  a  diameter  of  0.5  cm; 

ix)  saturation  spectral  and  line  radiances  were  calcu- 
lated from  expressions  given  by  Omenetto  and  Winefordner  (.40); 

x)  geometrical  factors  such  as  cell  sizes,  instrumental 
measurement  factors,  and  modulation  factors  were  chosen 
to  be  experimentally  realistic  for  atomic  fluorescence 
spectrometric  systems,  atomic  emission  spectrometric  systems 
(similar  to  commercial  flame  and  ICP  systems)  and  atomic 
absorption  spectrometric  systems  (with  a  narrow  line  source, 
similar  to  any  commercial  AA  spectrometer;  with  a  continuum 
source,  an  echelle  grating  spectrometer  was  assumed). 

Observations  and  Conclusions 

Several  conclusions  are  apparent  by  comparing  the  limits 
of  detection  (in  atoms/cm3  of  hot  gases)  and  the  S/N 
expressions  used  to  estimate  the  limits  of  detection. 

i)  Atomic  absorption  flame  spectrometric  detection 
limits  in  any  flame  or  plasma  are  in  the  range  of  106  -  109 
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atoms/cm  for  300  nm  resonance  lines,  and  1Q5  t   107  atoms/cm3 
for  600  nm  resonance  lines. 

ill  Using  a  line  source  in  atomic  absorption  spectrometry 
improves  detection  limits  anywhere  from  10X  to  50X  over  the 
use  of  a  continuum  source.  -• 

iii)  Atomic  emission  flame/plasma  spectrometric  detection 
limits  are  about  equal  to  those  of  atomic  absorption  spectrom- 
etry using  a  continuum  source  at  300  nm  (except  for  ICP)  or 
lower  wavelengths,  and  greater  than  50X  better  for  resonance 
lines  at  600  nm  or  higher;  ICP  typically  achieves  detection 
limits  in  the  range  of  103  -  105  atoms/cm3. 

iv)  Atomic  fluorescence  flame/plasma  spectrometer 
detection  limits  with  either  narrow  line  or  continuum 
conventional  sources  are  typically  107  -  109  atoms/cm3  for 
cases  where  "cell"  background  emission  shot  noise  is  dominant 
at  300  nm,  and  105  -  10s  at  600  nm.   If  the  "cell"  back- 
ground emission  can  be  made  negligible  and  only  dark  current 
shot  noise  remains,  then  detection  limits  with  either  line 
or  continuum  conventional  sources  will  be  102  -  105  atoms/cm3 
for  600  nm  resonance  lines,  or  104  -  106  atoms/cm3  for  300  nm 
resonance  lines;  the  ICP  and  the  H2/02/Ar  flame  are  much 
better  cell  environments  for  AF  because  of  their  much 
smaller  quenching  effects. 

v)  Atomic  fluorescence  flame/plasma  spectrometric 
detection  limits  with  laser  continuum  sources  and  saturation 
conditions   result  in  superior  values,  i.e.,  detection 
limits  in  the  101  -  102  atoms/cm3  range  for  the  N2-pumped 
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dye  laser,  in  the  1  -  10   atoms  /cm3  range  for  the  flashlamp 
pumped  dye  laser,  and  in  the  less  than  1  atom/cm3  range  for 
the  c¥  Ar*  dye  laser  assuming  detector  shot  noise  limitation 
and  approximately  1Q  ,  A  Or  and  101  atoms/cm3  for  background 
emission  shot  noise  limitation  for  the  three  laser  systems; 
background  noise  limited  detection  limits  at  600  nm  are 
always  better  than  or  equal  to  those  at  300  nm  because 
the  noise  level  is  generally  less  at  600  nm  and  signal  levels 
are  larger, 

vi)  Atomic  fluorescence  flame/plasma  spectrometric 
detection  limits  with  the  cw  Ar+  pumped  dye  laser,  under 
saturation  conditions,  are  better  than  with,  the  flashlamp 
pumped  dye  laser,  which  are  better  than  with  the  N2  -  pumped 
dye  laser  because  the  duty  factors  for  the  three  lasers  are 
0.5,  10~5  and  10"'  respectively.   However,  it  should  be 
noted  that  only  the  N2  -  pumped  laser  and  flashlamp  pumped 
dye  laser  beams  can  be  beam  expanded  to  approximately  0.5  cm 
diameter,  and  still  approach  saturation  (.see  typical  oper- 
ating conditions  in  Table  X  )  at  300  nm  or  even  at  600  nm. 

vii)  Atomic  fluorescence  saturation  limited  flame/plasma 
spectrometric  detection  limits  with  the  etalon  narrowed 
N2  -  pumped  dye  lasers  are  better  than  the  continuum  cases; 
it  is  also  clear  that  detection  limits  with  a  narrow  line 
dye  laser  decrease  linearly  as  the  square  root  of  the  line 
source  irradiance  increases.   Thus,  it  would  appear  that  if 
more  powerful  line  sources  were  utilized,  one  could  expect 
to  realize  a  further  enhancement  in  detection  limits.   In 
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addition,  in  the  case  of  resonance  excitation,  the  ratio  of 
fluorescence  to  fluorescence  scatter  will  he  considerably 
greater  for  a  narrow  line  laser  source  than  for  a  continuum 
laser  source  producing  the  same  fluorescence  signal  level. 

viii)  A.tomic  fluorescence  flame  spectrometry  using  a 
low  intensity  continuum  excitation  source  (flashlamp  pumped 
dye  laser),  in  an  C2H2/Air  flame  exhibits  a  detection  limit 
for  two  photon, (600  am)  absorption  via  a  real  level  one  half 
that  for  a  single  photon  (.300  nm)  absorption  under  similar 
conditions.   This  thereby  confirms  the  fact  that  the  two- 
photon  process  enhances  selectivity  only  at  the  expense  of 
sensitivity. 

ix)  For  the  case  of  two-photon  (600  nm)  absorption 
via  a  virtual  level  (no  saturation  possible)  the  gain  in 
selectivity  comes  at  a  much  higher  price  with  respect  to 

sensitivity.   In  this  case,  we  observe  a  limiting  observable 

6      7         3 
concentration  of  10   to  10   atoms/cm  ,  whereas  the  one  photon 

process  can  detect  from  10  to  104  atoms/cm3. 

x).  It  should  be  stressed  that  detection  limits  in 

3 

Tables  VITIandlX  are  in  atoms/cm  of  hot  gases  and  not  in 

solution  analyt'e  concentration  sprayed  (nebulized)  into  the 


cell.   If  the  process  of  conversion  of  analyte  into  species/ 
of  hot  gases  is  not  ef f icient,.itis  possible  for  a  technique 

such  as  laser  excited  AFS  in  a  flame  with  a  detection  limit 

3  3 

of  10   atoms/ cm   and  atomic  emission  ICP  with  a  detection 

limit  of  10   atoms/cm  to  invert  with  respect  to  detection 
limits  in  solution.   The  rather  impressive  detection  limits 


cm 


3 
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2  3 

(approximately  10   atoms/cm  )  by  laser  excited  AFS  have, 

however,  been  experimentally  approached  by  Fairbanks,  et  al. 
(_58)_,  by  Smith   et  al .  (591  and  by  Gelbwachs   et  al .  (60). 
For  the  theoretical  approach,  to  calculating  detection  limits 

with,  respect  to  analyte  concentration  in  solution  from 

3 

atoms/cm   in  hot  gases  see  Appendix  I. 

Analysis  of  Non-Photomultiplier  Techniques 
Atomic  Flame  Photoacoustic  Spectroscopy  (AFPS) 

The  final  two  techniques  which  we  wish  to  calculate 
theoretical  detection  limits  for  are.  photoacoustic  spectros- 
copy.  Although  these  two  techniques  are  not  highly  regarded 
as  acceptable  methods  of  practical  trace  metal  analysis,  a 
considerable  amount  of  research  is  being  focused  upon  them. 
Therefore,  it  seems  as  if,  besides  the  conventional  methods  of 
trace  metal  analysis  previously  discussed,  these  newer 
techniques  should  be  looked  at  with  respect  to  their  ultimate 
limits  of  detectability. 

To  begin  with,  we  will  discuss  the  analysis  of  gaseous 
species  by  AFPS.   Although.  AFPS  is  more  reknowned  for  the 
analysis  of  solids  (.61)  and  liquids  (62)  a  considerable  amount 
of  work  has  been  done  studying  the  analysis  of  gaseous  species 
(both  in  gas  cells  (.63,64,65)  and  flame  cells  (66)).   The 
approach  to  deriving  signal  and  noise  expressions  for  AFPS 
is  similar  to  that  taken  for  AAFS,  AEFS  and  AFFS.   However, 
due  to  fundamental  differences  in  the  modes  of  detection  of 
analyte  signal,  there  are  some  differences.   A  brief 
description  of  the  signal  producing  and  measurement  steps 
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may  prove  beneficial  to  understanding  the  final  equations. 

In  AFPS,  tiue  excitation  of  analyte  is,  analagous:  to  that 
of  AAFS  and  AFFS;  however  that  is  where  the  similarity  ends. 
In  AAFS,  we  are  interested  in  measuring  the  amount  of  light 
absorbed  by  an  analyte,  in  AFFS  our  attention  is  focused  on  the 
amount  of  light  fluoresced  by  the  analyte  and, in  AFPS ,  we 
focus  our  attention  on  the  translat ional  energy  produced  by 
collisional  de-excitation  of  the  analyte  species.   The  first 
two  techniques  establish,  quantitative  relationships  between 
analyte  concentration  and  light  energy,  wheras  AFPS  relates 
translat  ional  energy  (pressure),  to  concentration.   The  detection 
device  most  commonly  used  in  AAFS,  AEFS  and  AFFS  is  the 
photomultiplier  tube,  whereas AFPS  utilizes  a  suitable  pressure 
transducer,  namely  a  microphone.  * 

We  are  interested  in  measuring  the  signal  generated  by 
laser  excitation  of  atoms  introduced  into  a  typical 
air/acetylene  flame.   This  will  be  accomplished  by  monitoring 
the  pressure  wave  resulting  from  the  periodic  analyte 
excitation  with  a  condenser  microphone.   The  analysis  is  based 
on  the  following  assumptions: 

i)_  all  atoms  that  absorb  the  incident  radiation  will 
either  return  to  their  original  ground  state  via  radiational 
deactivation  (fluorescence)  or  radiationless  deactivation 
(.quenching  collisions).; 

ii)_  those  atoms  that  return  to  the  ground  state  via 
quenching  collisions  will  do  so  by  transferring  their  energy 
into  translational  energy  of  the  flame  gas  molecules  (.67); 

iii)_  this  periodic  energy  conversion  process  (whose 
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frequency  is  determined  by  the  frequency  of  the  chopped  light 
source!  will  in  turn  produce  a  pressure  waye  whose  magnitude 
will  remain  constant  until  detected  by  the  microphone  which 
is  located  approximately  4.  cm  away  from  the  flame  border; 

iv)_  the  absolute  limiting  noise  source  is  taken  to  be 
the  equivalent  noise  pressure  impinging  upon  the  diaphragm 
of  the  microphone  which  is  due  to  the  Brown i an  motion  of 
the  local  air  molecules. 

Although  this  final  assumption  does  not  fairly  represent 
the  limiting  experimental  noise  sources,  e.g.,  environmental 
noise  pickup,  chopper  vibration  pickup,  flame  background 
noise,  flame  temperature  fluctuations , etc. ,  it  is  assumed 
that  with  sufficient  care  and  effort  these  experimental 
noises  can,  in  principle,  be  eliminated  leaving  the  limiting 
noise  source  to  be  that  due  to  Brownian  motion. 

As  mentioned  for  the  cases  of  atomic  spectroscopy 
using  a  photomultiplier  detection  schene ,  the  following 
results  will  reflect  the  absolute  best  conditions  possible 
for  a  typical  photoacoustic  experiment  and  will  express  the 
lowest  possible  detection  limit  the  analyst  can  hope  to 
achieve.   The  following  is  the  derived  expression  for  the 
signal  produced  (_in  units  of  energy  per  unit  volume), 

4 
s  =  EcaoX  h        fa  ^-Y>k21  nL  S  Ql  )  £  x. 
I  8ttc  *1  '         T    47T     X 

where  all  the  terms  have  been  previously  defined  except  for, 
k21   =  collisiona-l  rate  constant  for  radiat ionless 
deactivation,  s~ 


87 

8  h  fraction  of  collisional  energy  that  is.  transferred 
to  transnational'  energy  of  th.e  gas  molecules;, 
dimensionless. 

ft  s  fraction  of  4tt  sr  collected  by  the  microphone, 

477 

dimensionless 
The  limiting  noise  expression,  as  previously  mentioned, 
will  be  derived  from  the  Brownian  motion  of  molecules- 
impinging  upon  the  diaphragm  of  the  microphone  and  is 
expressed  as  follows, 

N  =  /  8k  T.  m  Af  \i 

m    * 
where, 

k  =  Boltzmann's  constant,  erg  K 

T  =  temperature,  K 

m  -   mass  per  unit  area  of  the  diaphragm,  g  cm-2 

Af  =  frequency  interval  about  the  measured  frequency,  Hz 

M  =  area  of  diaphragm,  cm 

Tm  =  damPin&  time  of  the  diaphragm,  s 

In  order  to  evaluate  the  S/N  expression,  numerical 
values  for  the'  variables  will  be  the  same  as  those  in 
Table  X  except  where  specified  otherwise.   The  calculations 
were  performed  for  atoms  whose  wavelength  of  excitation 
is  centered  at  3Q0  nm  and  at  600  nm.   We  will  look  at  a 
"typical"  good  absorber  and  weak  fluorescer.   The  light 
source  is  the  cw  Ar+  pumped  dye  laser.   The  microphone  is 
represented  by  a  B  &  K  4144  type  condensor  microphone.   A 
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10.  cm  pathlength  was  chosen  for  the  flame  cell,  while 

-3 

a  total  of  5  X  10.  *  of  the  total  amount  of  sound  produced 

is  collected  when  considering  the  geometry  of  the  microphone. 

The  results  of  the  calculations  are  given  in  Table  XI. 
With  respect  to  typical  limits-  of  detection  reported  in  the 
literature  ranging  from  1  X  107  to  1  X  1010  species/cm3 
(.68,69),  the  results  appear  to  he  very  reasonable. 

Atomic  Flame  Optogalvanic  Spectroscopy ( AFOS) 

Atomic  flame  optogalvanic  spectroscopy  has  been  shown 
to  be  another  valuable  method  for  measuring  trace  metals  in 
flames  (.70,71,72).   The  fundamental  principle  of  AFOS  is 
the  ionization  of  the  analyte  species  and  monitoring  the 
current  produced  when  the  cell  is  in  the  midst  of  an 
applied  voltage.   In  the  derivation  of  the  expression  for  the 
signal  produced  by  AFOS,  the  experimental  system  used  as  a 
model  was  that  of  Travis  et  al .  (73).   They  utilized  a 
flame  cell  (air/acetylene)  10  cm  long  equipped  with  two 
parallel  molybdenum  plates  as  the  ionization  current 
measurement  device.   The  source  was  a  flashlamp  pumped 
dye  laser  with- a  pulse  width,  of  1  X  10~6  s.   The  limiting 
noise  source  was  taken  to  be  shot  noise  due  to  the  inherent 
electron  population  of  the  flame.   Four  cases  are  to  be 
examined,  all  of  which,  are  calculated  for  a  typical  atom 
whose  absorption  line  is.  centered  around  300  nm.   The  first 
two  are  those  in  which,  the  atom  has  a  high  ionization  rate 
constant  (74)  and  high  and  low  background  noise  levels.   The 
last  two  cases  are  for  a  low  ionization  rate  constant  and  the 
same  varying  background  levels. 
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Tattle  XI 


Calculated  Detection  Limits  for 
Atomic  Flame  Photoacoustic  Spectroscopy 


Wavelength         Concentration 

(nra)  3 

(atoms /cm  ) 

300  2.0  X  109 


600  1.3  X  107 


90 
Table  XII 


Theoretical  Limits  of  Detection  for 
Atomic  Flame  Optogalvanic  Spectroscopy 


3 

Absorption  Process         Concentration  (atoms/cm  ) 


One  photon 
:2i 


low  background    high  background 


k„.  -  108  s  1  5  X  106  5  X   107 


k2i  =  104  s_1  5  X   1010  5  X  1011 


Two  photon 


19-1  3  A. 

virtual  Qs.  ±   =  10   s   )  2  X  10  2  X  10* 

19   -1  -=>  -4 

real  (k   =  lCT"5  s  x)  3  X  10  °  3  X  10  * 

O  1 
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The  intensity  expression  is  given  as  follows  (JLn  eyents), 
S  =  VpCA^o1'^^^  Y  Attk21  Ti 

8TTC  g. 

where  all  the  variables  have  heen  previously  defined  except 
for  the  following, 

I     =  length  of  collector  plates,  cm 

M  -  height  of  collector  plates,  cm 
k_.  =  ionization  rate  constant,  s 

The  background  limiting  shot  noise  Cin  events),  due  to 
the  electron  population  of  the  flame  is  calculated  via  the 
following  expression, 

N  "  ^2DD  *B  Ti^ 


where , 

R„  =  average  background  count  rate,  events/s 

D 

If  we  assume  a  background  current  of  1  X  10  A  for  a 

-6 
low  background  flame  and  1  X  10  A  for  a  high  background 

flame  C41),  a  10  s  integration  time  and  that  the  high 

ionization  rate  constant  is  1  X  10   s~   and  the  low  ioni- 

4   -1 
zation  rate  constant  is  1  X  10   s   ,  we  are  now  able  to 

calculate  the  theoretical  limits  of  detection  for  AFOS . 

The  results  can  be  found  in  Table  XII.   These  results 

agree  extremely  well  with  those  of  Travis  et  al.  (7  3)  who 

estimated  that  their  experimentally  derived  detection  limits 

ranged  from  1  X  10   atoms/cm  to  1  X  10    atoms/cm  . 
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In  an  analagous  fashion  to  atomic  fluorescence  flame 
spectrometry,  AFOS  has  been  shown  to  gain  certain  advantages 
through  the  utilization  of  two-photon  absorption  processes. 
As  in  both,  cases,  the  two  photon  absorption  process  enhances 
the  attractiveness  of  the  techniques  by  enhancing  their 
selectivity.   However,  due  to  the  relatively  weak  absorptivities 
in  the  two  photon  process  (via  a  virtual  level)  this 
enhancement  cones  only  at  the  expense  of  sensitivity. 
Nevertheless,  we  shall  consider  the  cases  of  two  photon  (600 
nm)  absorption  via  a  virtual  level  and  that  via  a  real 
level  for  AFOS  (Tigure  II).   For  the  case  of  the  two  photon 
absorption  via  a  virtual  level,  the  S/N  expression,  assuming 
the  same  limiting  background  noise  source,  becomes , 

2  L 

2  photon,  virtual   En    AxT  <5  Y  n„  A  I   k„ .  t. 

S/N)  =   LP    L  l  __ 


lim,  BS 


i 


(2DD  RB  r.r 


For  the  case  of  a  real  intermediate  level,  the  S/N  expression 

is  as  follows, 

L 
2  photon,  real    Dq  nT  A  £  k  .  x. 
S/N)  -    b  ~4  31   x 


lim,  BS 


i 


(2DD  RB  t.)2 


where  all  terms  in  the  preceding  two  equations  have  been 
defined  except , 

Att   =  pulse  width  of  the  laser,  s 

k„.   =  ionization  rate  constant  from  third  level,  s~ 
The  term  k3i  is  derived  from  k„ .  by  the  Boltzmann  relation- 
ship, 


93 


ionization  limit  45,000  cm 


-1 


-1 


k3i 
3                         33334  cm 

2 

k2i 

16667  cm 

-1 


Figure  II 


0  cm 


-1 


Energy  Level  Diagram  for  a  Hypothetical  Atom 
Undergoing  Photoionization 
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dl    Z1  expC-AE2./kT) 
where  £E„ -  is  the  energy  difference  Between  level  three  and 
the  ionization  threshold  (Tigure  II)  and  AE„ .  is  the  energy 
difference  Between  level  two  and  the  ionization  threshold. 

The  results  of  the  calculations  are  to  be  found  in 
Table  XII.   With  respect  to  the  conventional  one  photon 
process, a  two- photon  absorption  via  a  virtual  level  exhibits 
a  detection  limit  three  orders  of  magnitude  better.   However, 
the  ultimate  in  sensitivity  is  seen  for  the  two  photon 
absorption  via  real  levels.   This  is  mainly  due  to  the 
assumption  of  saturation  of  each  transition  which  is 
theoretically  appealing ,  although  experimentally  difficult 
to  obtain. 


APPENDIX  A 

DETECTION  LIMITS  BY  ATOMIC  SPECTROMETRY  METHODS 
FOR  ANALYTES  fg  REAL  SAMPLES 

As  previously  discussed,  in  analytical  atomic  spectros- 
copy, it  is  necessary  to  convert  the  analyte  in  a  sample 
into  gaseous  atoms  (or  ions)  in  order  for  atomic  processes 

to  be  observed.   One  is  therefore  interested  in  the  conversion 

3 

factor,  9,  from  atoms/cm  of  hot  gases  to  ug/ml  of  sample 

solution  assuming  the  sample  is  nebulized  to  produce  the 
atomic  vapor.   The  nehulization  process  is  the  most  common 
means  of  sample  introduction  although  electrothermal  atomizers 
have  widespread  use  as  well.  However,  for  reasons  previously 
mentioned,  they  will  not  be  discussed  here. 

The  conversion  factor,  Qr    is  given  by  the  following 
equation  (41) , 

c    10.  Q  e  M      Ug/ml 
0  =   L  = 


nT   N.  F  s  |3   B  ■   atoms/cm3 

where  the  terms  are  defined  as  following; 

c  =  solution  concentration  of  analyte,  yg/ml 

L 
n  3 

T   =  limiting  gaseous  concentration  of  analyte,  atoms/cm 

MA  =  atomic  weight  of  analyte,  g/mole 

F  =  transport  rate  of  sample  solution  as  determined 

3 
by  the  nebulizer,  cm  /s 
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e  -  efficiency  of  nebulization  and  desolvation,  dimen- 

sionless 

g  =  efficiency  of  vaporization,  dimensionless 

3   =  efficiency  of  atomization,  dimensionless 
a 

Q  =  flow  rate  of  unburnt  gases  (at  room  temperature) 

3 
into  the  flame,  cm  /s 

ef  =  flame  gas  expansion  factor,  dimensionless. 

Therefore,  to  convert  the  estimated  detection  limits  present 

3 
throughout  this  work  (_from  atoms/ cm  to  ug/ml)  one  just  needs 

to  multiply  the  tabular  values  by  0. 

A  brief  discussion  concerning  the  variables  involved 

may  serve  to  exemplify  their  relationships  to  the  0  factor. 

For  instance,  e  can  vary  from  approximately  0.01  to  0.1 

depending  upon  the  nebulizer  chamber  design  (.7'5)  but  is 

nearly  independant  of  the  flame/plasma  being  used  for 

analysis.   On  the  other  hand,  %      depends  critically  (especially 

for  cooler  flames  such  as,  C„H_/Air  or  H_/0~/Ar)  upon  the 

sample  matrix,  i.e.,  the  chemical  components  and  their 

relative  concentrations  to  the  analyte  concentration,  and 

upon  the  analyte  chemical  and  physical  properties.   For 

hotter,  more  reducing  combustion  flames,  as  the  C^Hp/NpO 

flame,  and  especially  for  the  hot  inert,  ICP  discharge,  these 

problems  are  much  less  severe.   Finally,  the  parameter  6 

depends  critically  upon  the  thermodynamic  characteristics 

of  the  analyte,  the  concentration  of  radicals,  as  0,  OH, 

H,  CN,  CH  etc.,  in  the  flame  gases,  and  on  flame  temperature. 

Parsons  et  al .  (48)  have  tabulated  0   values  measured  by  a 

a 
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number  of  researchers  for  a  number  of  elements  in  several 

flames. 

In  order  to  estimate  0  values  and  the  corresponding 

solution  detection  limits  from  flame  detection  limits,  a  few 

assumptions  must  be  made.   If  we  assume  an  CnH  /Air  flame 

2  2 

produced  by  2  1/min  of  C-H-  and  10  1/min  of  air,  then  Q 

3 

is  200  cm  /s  of  flame  gases  at  2500  K  Of  =  8.5).   Assuming. 

3 

a  transport  rate  of  3  cm  /min  and  an  efficiency,  z,    of  0.10, 

3 

then  F(ef)  is  0.005  cm  /s.   Therefore, 

6.  X  10~13M. 

0  =  a  . 

3  6 
v  a 

For  Mg  CM  =  23  amu) ,  8   B   M  and  so, 

A  y    St 

0(Mg)  =  2  X  10"11 

6  3 

Thus  a  detection  limit  of  10   atoms. /cm   (as  in  AA)  would 

become  0.02  ppb  (ng/ml)  in  solution.   If  these  calculated 
results  are  compared  to  experimentally  derived  limits  of 
detection  tabulated  in  Table  XIII  good  agreement  is  observed. 
For  instance,  the  calculated  limit  of  detection  of  0.02  ppb 
for  Mg  by  atomic  absorption  spectrometry  with  a  line  source, 
is  only  one  order  of  magnitude  below  that  which  is  experiment- 
ally achievable.   Furthermore,  if  we  assume  0  to  be  10~ 
for  any  typical  atom, then  the  theoretically  calculated 

atomic  emission  flame  spectrometry  detection  limits  of 

—2       —8 
10  '"'  to  10  '  ppm  favorably  compare  with  the  experimentally 

1      -4 
obtainable  concentrations  of  10  to  10   ppm. 
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ABBENBIX  B 


DEFINITION  OF  FLICKER  FACTOR 


The  flicker  factor,  g,  can  easily  be  defined  if  one 
considers  the  spectral  noise  power  for  flicker  noise 
which  is  giyen  in  terms  of  current  fluctuations  by, 

SflC«  =^~  -   k2  ^  *2 

f  f 

where, 

f  =  frequency,  Hz 
K  =  a  constant  describing  the  low  frequency  stability 
of  the  noise  source,  dimensionless 
i  =  current  associated  with  the  noise  source,  A 
R  =  measured  rate  associated  with  noise  source, 

counts /s 
e  =  elementary  charge  of  electron,  coulomb 
Having  determined  K  one  can  define  the  flicker  factor 
£,  for  integrator  measurements  as, 
£.  ■  CL83  K. 
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